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OlAPTER 1 


INTRQIXJCriQN 

Since helicopters anci other VTX aircraft are designed to take-off and 
land vertically, to hover for rescue attempts and other purposes, and also 
to be able to fly close to the ground for a long time, study- of their per- 
formance in ground effect is very important. 

Calculation of stability and control derivatives of helicopters and 
aircraft requires knowledge of the behavior of aeroi^namic forces and mo- 
ments for different flight conditions. Since the wake of fixed-wing 
aircraft is left behind in forward flight, there is no interaction between 
the newly generated vortices and old vortices. This simplifies the wake 
stu<^ of fixed-wing aircraft compared with helicopters. There are three 
major differences associated with helicopter wakes. 

First, the wake of a helicopter does not move away from the rotor as 
it does in conventional airplanes. In slow forv/ard flight, the vortices 
shed from the blade tips move downward belcw the rotor as a helix. 
Secondly, the problem becomes more complicated as the helicopter approaches 
the ground. In the presence of the ground the wake contacts the ground, 
rebounds and in certain flight conditions is drawn through the rotor again. 
Thirdly, the wake can roll up in front and around a helicopter like a 
horseshoe. This vortex (which is close to the ground) is called a "ground 
vortex". As speed increases a helicopter overruns the ground vortex ahead 
of it causing a transient disturbance well known to pilots. 

All these problems complicate the stucfy of helicopter aerodynamic 
forces and moments. Most of the investigations of the helicopter wake in 
ground effect have been experimental and qualitative, although there has 
been seme theoretical study- of hover in ground effect by modeling the wake 
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as flat vort'ix rings* So far, there has not been a sui-aessful stud^ of the 
wake in ground effect which covers all flight conditions because of dif- 
ficulties irentioned previously. The current stu<V has been made to address 
these problems. 

Precise mathenatical description of the flow field for a rotor is made 
difficult by the interdependence of the velocity and the wake position. 
That is, the boundary condition is known but the location of the boundary is 
not known. This situation is characteristic of problems that can be 
generally classified as "free-boundary problems". The method of solution 
that can scmetimes be used on this class of problens is to guess the posi- 
tion of the boundary, compute the solution, and determine if the computed 
solution is consistent with the assumed boundary location. This approach is 
based on the use of space coordinates as independent variables; computation 
j,s cast as some sort of feedback scheme in which the differences between the 
assumed boundary location and the resulting solution are used to determine a 

new estimated boundary position. Advantage is taken of the azLmuthal peri- 

0 

odicity of the rotor position and independent variable is chosen to be time 
(or azimuth angle of a reference blade). The point of view, then, is that 
new elanents of wake vorticity are generated as the blades rotate and 
translate. These elements are onvected with velocities determined by 
their self-induced velocities as well as velocities induced by existing wake 
structure, the bound vorticity of the blades, and the iiiage vortices intro- 
duced to satisfy the gicuiid-plane boundary condition. It is assumed that 
eventually a stabilized periodic wake array solution can be obtained, since 
the condition of a fixed periodic blade loading was imposed. 


The continuous distorted helix was selected as the model for the cur- 
,cant effort. The ground boundary condition was enforced by influence of an 
image wake. 

In the present study it has been found that the distorted helix model, 
in the absence of ground yields average velocities that agree well with 
measurenents made in and about a helicopter rotor wake (reference [2]). 
There is a good qualitative check on the wake of a helicopter in ground ef- 
fect, but no check 'lould be made quantitatively because these measurfments 
are not available. 

The digital computer program is based on the assumption that the rotor 
is in stea<^ level flight (or hovering) with a specified tip-path-plane 
angle. Shaft rotational speed, rotor force, initial vortex core radius, 
number of blades and the ratio of the rOtor radius to height above the 
ground are also inputs. Information relative to the computing program is 
given in Appendix B. 

It is the purpose of this stucfy to provide the numerical methods from 
which to calculate aerodynamics (stability derivatives) necessary for the 
modeling of helicopter needed for design of stability and control 
^stans. 


CHAPTER 2 


IDEALIZED MODEL 

Conputations of aerot^namic forces, moments, and stability derivative 
of helicopter directly depend upon ir»3uced velocity at the rotor disc. As 
induced velocity is a consequence of a wake structure, it is necessary to 
have a wake structure which gives relatively accurate velocity on the rotor 
disc. 

In addition it is important to make the model as simple as possible to 
reduce expensive computation time. In the following sections a wake model 
will be presented and studied. Section 2-"l describes and compares the 
prescribed and free wal;e analysis. Sections 2-2 through 2-7 explain the 
free wake structure and discuss the reason for keeping important parts of 
the wake arid ignoring unimportant parts of the wake. In section 2-8 the 
basic formula to calculate induced velocity of a rectilinear vortex segment 
at a point in the space will be formulated. Section 2-9 describes the 
structure of the computational model, and in section 2-10 the induced 

t 

velocity formula will be modified for special cases. Section 2-11 considers 
the inboT:rd wakes and firoally discusses the reasons for considering only one 
of tlie inner wakes. 

2-1 FREE WAKE ANALYSIS AND PRESCRIBED WAKE ANALYSIS. 

Wake analysis of helicopters has been the topic of many researchers 
for many years. It is one of the most important aspects of helicopter 
study, because the majority of the characteristics of a helicopter depends 
on its wake structure. 

.Among the many metliods of analyzing helicopter wakes, two methods cur- 
rently are employed by the majority of helicopter researchers; prescribed 
wake analysis and free wake analysis. 
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In the former method, for each flight condition the location and in- 
t<iU)sity of the vortices of a wake are measured. Having these locations and 
intensities, aeroc^namic forces and moments are computed. As previously 
mentioned, for each flight condition a number of experiments for different 
points in the space must be carried out. An increase in accuracy can be ob- 
tained by increasing the number of points on the wake to be examined. 

In the free wake analysis, an initial wake with seme initial 
properties is assumed. Ihen it is possible to compute the induced velocity 
of any point in the space, including the points on the wake. The new loca- 
tion of the wake solution can be obtained after an increment of time by 
using velocities of the points on the wake with simple forward integration. 
Also the properties of the wake can be mcx3ified to be consistent with tiie 
new wake. This procedure can be repeated until the wake location stabilizes 
and reaches its periodic steac^ state. 

Similar to the prescribed wake analysis, free wake analysis accuracy 
c:an be increased by increasing the number of the segments. Because the com- 
putation time will grow exponentially with increase in number of segments, 
there is a practical limit on the number of segments. 

To study the wake in ground effect, the detailed analysis of the 
prescribed wake is prohibitively expensive in terms of manpower required and 
difficulties in measuring the location and circulation of the wake segments 
close to the ground, especially in presence of a ground vortex. Free wake 
analysis is considered to be a better alternate for the wake analysis, be- 
cause it can cover all the points in the wake including the points close to 
the ground in almost all flight conditions, and it is more efficient and 
less expensive. For these reasons the free wake method was chosen to stu(^ 
the ground effect. 
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2-2 WAKE SfIRUCrURE. 

In a given flight condition, the wake of a helicopter contains 
several different vortices. Considerir*5| potential flow around an airfoil, 
each blade nay be replaced by a vortex line having the same lift and ap- 
proximately the same flow field. This vortex line is called a bound vortex. 
A complete model of a bound vortex for each blade consists of radial varia- 
tion of circulation, as well as tangential variation with azimuth angle. 

Due to steep change in circulation at the tip of a blade, tip vortices 
are generated; and because of the variation of vorticity along each blade, 
inboard vortex filaments parallel to tip vortices are released. Blade pitch 
angle variation or velocity change of a blade results in vortex shedding by 
the rotating rotor blades (figures 2-1 and 2-2) . 

Near the ground, a mirror image of the whole wake must be added to 
satisfy the boundary conditions on the ground. 

Consideration of the ground effect will result in large increases in 
computation time for the following reasons: 

i) Computation of the induced velocity of the mirror image wake doubles 
the computation time. 

ii) The number of points for wake analysis in-ground-effect should be 
more than out of ground effect wake study. Hie reason for this increase is 
that during flight out of the ground effect, as time passes the vortices 
move away from the rotor and their effects can be neglected; whereas in 
ground effect, old vortices may hit the ground, bounce back and in some 
flight conditions interact with the wake and cause major variations in 
forces on the rotor. Therefore it is very important to keep the effect of 
the old wake which can have a significant effect on the velocity distribu- 
tion of the rotor. 
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Figure 2-1* Schematic of Rc:cr '.•.'ake Structure. 







Figure 2-2. Segmented Discrete Vlortex Rapresentatioa of the Wahe 



A detailed analysis of the whole wake is practically iitipossible espe- 
cially in the presence of the ground and ground vortex. Therefore, it is 
important to make the model as simple and as accurate as possible utilizing 
the parts which describe the wake efficiently for velocity computations, 

2-3 TIP AND ROOT VORTICES 

As was mentioned in the previous section, a quick change in circula- 
tion close to the tip of the blade causes a steep continuous sheet of 
vorticity close to the tip to be released (see figure 2-1 to 2-3 for con- 
tinuous sheet vortices and discrete filaments) . Flow visualization studies 
indicate that this sheet of vorticity rolls up within a few chord lengtiis of 
the blade and forms a single concentrated vortex line. To avoid complexity 
it will be assumed that tip vortices are fully rolled up from beginning as 
they are released. Tip vortices are the most important part of the wake, as 

p%s:r,f a considerable amount of energy. Ihe velocity distribution of 
the rotor disc greatly depends on the locations and strengths of the tip 
vortices. Therefore the primary goal is to compute the locations and 
strengths of the tip vortices. 

Similar but weaker vortices are created near the blade root, however 
analytical and experimental studies have shown that they rapidly dissipate. 
Even if root vortices are considered, their contribution to the induced 
velocity is negligible. For this reason, as well as the savings in computa- 
tion time, they will not be considered. 

2-4 BOUND VORTICES 

Among many models available for bound vortices, the simplest one has 
been chosen to replace the blades and their images. Each blade is replaced 
by a single radial vortex line with constant vorticity along its length. 
The strength of these vortices is computed by assuming the total load on all 


tne blades afproxinately ijquals the helicopter weight. The reason for chos- 
ing such a simple mo^el is to avoid large increase in computation time. 

2-5 imOARD VORTICES 

As previously mentioned, the accuracy will increase with utilization 
of more realistic models. Tlie mathematical model will be nvore realistic if 
large number of inboard vortex filaments are included. Consideration of 
even a few inner vortices will result in an increase in computation time by 
a large factor. Therefore it is desired to include as few inner wake fila- 
ments as possible. It will be shown that if only one inside vortex filament 
at r/RrO.7 is included, the accuracy of the velocity close to rotor plane 
will iirprove considerably. Also computation time is only increased 
slightly. 

2-6 SHED VORTICES 

In forward flight, the pitch angle of each blade will vary peri- 
odically (once per revolution) . This variation will result in variation of 
circulation of each blade with azimuth angle of the rotor. Consequently, 
periodic variation in blade circulation results in release of a continuous 
vortex sheet parallel to the blade which may be modeled as vortex filaments. 

Fortunately, during flight conditions in which ground effect is important, 
the vortex filaments are not very strong and they can be neglected. 

2-7 MIRROR IMAGE VORTICES 

To satisfy boundary conditions on the ground, it can be assumed that 
corresponding to each blade bound vortex and each wake segment, there is an 
image vortex with opposite circulation belcw the ground at the same 
distance. The induced velocity of the image segments at points close to the 
rotor disc is snail and for the points close to the ground, the induced 




velocity is large. This effect forces the wake to satisfy the boundary con- 
dition on the ground and consequently forms the v;ake close to the ground. 
In sane flight conditions the wake close to the ground rolls up, passes 
through the rotor and forms a new configuration for the wake (ground 
vortex) . Iherefore, it is necessary to take into account the image wake as 
an important part of the whole wake. 

2-8 EX)RMULATION OP IHE ^DDEL 

Formulation of the flow corresponding to the simplified model can be 

accomplished as follows: A coordinate system fixed in the tip path plane, 

a plane which passes through tip of the blades, is introduced (figure 2-4) . 

The vector V in Figure 2-4 is the free-stream velocity (the negative 
f 

of the translational velocity of the aircraft) , inclined at an angle a 

T 

to the X-axis and parallel to the X-Z plane. Ihe angle ip denotes the 

azimuthal positioning of a given point with respect to the origin. 

The air velocity V at a given point located by the vector r riiay be 

P P 

expressed in the form 

nr )- -1-f ( 2 - 1 ) 

P 

ly 

where r =r - r and r- ( r ) is the circulation about the vortex element at 
1 P 

r. The line integral is to be taken over all vortices in the flew and the 
image systan whose paths are collectively denoted by C . It is necessary to 

modify this expression when r =r and this point is .scussed later. The 

P 
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circulation/ r t about the blade vortices is estimated directly in terms of 
flight parameters. The circulation about a wake vortex at a given point is 
simply that value assigned to the blade vortex when it was generating that 
segment of the wake. 

The only information lacking for the complete specification of the 
flew at a given instant/ then/ is the location of the wake vortices at that 
instant. The position of a given point on a wake vortex located by the vec- 
tor r is the time integral of the velocity experienced by that fluid 
particle 




r(t)=r(tj+# V(r(r)) dr 

't. 


( 2 - 2 ) 


ThuS/ even after being simplified/ the flew can only be obtained as 
the solution of the nonlinear integral formed by the substitution of 
Equation (2-1) into Equation (2-2). A direct analytical solution is not 
feasible/ but the problem is amenable to solution by numerical methods using 
a high speed computer. The manner in which the formulations of Equations 
(2-1) and (2-2) were in^sleraented for digital computation will be discussed 
in the following sections. 

2-9 DESCRIETION OP CX5MTUTATI0N MODEL 

For the purpose of numerical analysis/ the wakes produced by each 
blade are divided into small segments (figure 2-5) . These segments are 
chosen to be sufficiently short so that/ they can be considered rectilinear 
vortices having constant circulation along their lengths for purpose of com- 
putation of the induced velocities. 
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Each segment is defined by its two ends. Segment S is between point 

i 

i and point i+1. The wake configuration at any instant is defined by the 
location of these end points. 

Computation is initiated ky specifying that each of the wake vortices 
lies on a prescribed curve. The curve is chosen to be a single vertical 
helix, or a previous solution for. a flight condition close to the required 
one. The numerical version of equations (2-1) and (2-2) are then performed 
by first summing the velocity contribution of all the vortex elements in the 
flew at each reference point (Equation 2-1 ) , and then using these velocities 
to compute the new location for each point for a time interval d t 
(Equation 2-2). The time interval d t is chosen to correspond to a small 
finite change in the azimuth position ^ of the blades. 


Axlr^Q dt 


(2-3) 


where Q is rotor angular speed. Once the new coordinates of each reference 
point is computed, the azimuth of the blade vortices is increased by^ and 
the velocity computation is performed again. As the blade vortices are 
repositioned, a new wake vortex element is added to the flew at the tip of 
each blade vortex, the added vortex having a length of approximately R d , 
where R is rotor radius. A corresponding element is dropped from the com- 
putations at the downstream end of each wake vortex to maintain a wake of 
constant size. The computations are continued in this manner for a 
specified time and the results are inspected. If a nearly periodic solution 
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is established in the space volume of interest, the calculation is 
terminated. 

The total number of the wake vortices taken into account and the mag- 
nitude ot A \J/ determine the accuracy of the flow representation at a given 
point. It is believed that, for a two-blade rotor, a value for d of 
thirty degrees is sufficiently small to furnish an acceptable estimate of 
the time variations of the flow 'consistent with the other approximations 
introduced. The number of wake elements to be considered depends on the 
region of interest, forward speed and height of the rotor above the ground. 
If the free stream does not clear the wake under the rotor, the number of 
wake elements must be sufficiently large to include all the wake elements 
close to the rotor. This phenomenon greatly depends on forward speed and 
height of the rotor above the ground. 

2-10 VELOCITY INDUCED AT POINT P BY THE AND THE BLADES. 

The velocity induced at an arbitrary point P by the vortices repre- 
senting the wake and the blades is simply the sum of the effects of an array 
of rectilinear vortex segments. If V denotes the velocity induced at P by 

the elanents between points P and P , it is found from equation (2-4) (Ref. 

1 2 

II] Page 152) that 


V = r (Cosd -Cos Q )/(4 7T h) 
1 2 


(2-4) 




where r is t±e strength of the element and 0 , $ and h are defined in 

1 2 

figure 2-6. The velocity is directed nocnnal to the plane containing P , 

1 

P and P. 

2 

As the field point P is snade to approach any point on the line joining 

P and P the induced velocity increases without limit, because h tends to 
1 2 

zero (Equation 2-4); the velocity becomes indeterminate for h=0. Because 
the velocity of air can not reach infinity, another model is employed for 
anall h. 

Experimental studies [2] of the structure of trailing vortices shcv; 
that for sufficiently snail h, the flow rotates as a rigid body. The region 
where vortex filanents have rigid rotation is called the core. Among many 
models suggested for r of the core, Scully's model [3] is believed to give 
the best results. 



(Vo) ^ 

I (h/pj2 


(2-5) 


Here a is defined as the radius of the core. The Scully model approaches a 
potential flew just a few core radii away from the filament. Because of 
favorable comparison with experimental data and snoothness for small h [4], 
the Scully model was used in all calculations. Figure 2-7 compares the nor- 
malized vorticity and normalized velocity of the Scully model with wind 
tunnel experimental data [2]. One of the advantages of this model over 
other models is that there is not a sharp change in velocity profile at the 
boundary of the core. This smooth behavior helps to avoid large changes in 
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computation of the induced velocity at a point close to a vortex filament 
which may result in numerical instability. 

In previous formulations, point P was not assumed to be one of the end 
points of vortex segments. If point P is one of the end points, ti^en it 
should be assumed to be a point on an arc segment containing the two seg- 
ments which P is one of their ends and undergoes self-induced velocity. 

’I » 

The self-induced velocity of the two segments is formulated in appedix 

A, The velocity of the point P is considered to be the same as self- 

2 

induced velocity of an arc consisting half of S and half of S . 

1 2 

V «rrin((Wa )*Tan((^ /4) )-1.]/(87TR)-<* 

^ 1 1 

r[Ln((av/a )^Tan( 0 /4))-l.l/(8n'R) (2-6) 

2 2 2 

where <f> and ^ are defined in figure 2-8. The self-induced velocity is 
1 2 

directed normal to the plane of the arc of the two segments. The ap- 
proj;lmate core radius of a given elanent may be assigned on a rotational 
basis using energy ^considerations. 

2-11 INNEK WAKE 

For computation of the stability derivatives the induced velocity on 
the rotor disc should be known more accurately than at the points far from 
it. To improve on simple model of the blade with constant vorticity along 
its length, the blades are discretize into tiny segments and all small in- 
board wake filaments parallel to the tip vortices are considered. Since 
computation time limits the total nuriiber of the segments only one inside 
wake at r=0.7R with =0.5 was assumed. Since the slope of the circula- 
tion along the blade is steeper around r=0.7R (Figure 2-9) consideration of 
an inside wake at this point is a better representation of the wake. It is 
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expected that this mcdel will have better results than any other single in- 
side wake. 









CHAETER 3 


KUt>ERICAL FQRMUTATTON 

The models of the rotor, the wake and image v^ake described in the pre- 
vious chapter were formulated as oountinuous functions of time. A digital 
computer can not integrate continous function exactly, therefore step-wise 
and interpolative approximations have to be made. 

A rectangular integration schgite is used in performing integrations in 
time. lhat is, when integrating velocity to compute displacanent, the 
velocity is assumed to rauain constant over the interval of time correspond- 
ing to a snail finite change in the azimuth position of tlie blades. 

Special integrations over the wake vortices are performed by assuming 
that these vortices are made up of snail rectilinear vortex segments whose 
circulation is constant form one point to the next. The position of the 
wake is then defined by the location of the end points of these segments. 
Consistent with the approximation made in the time integration, the initial 
length of each wake segments the length of the arc swept out ty the blade 
tip over the interval used for time integration. Self-induced effects at a 
given wake point are computed by taking, as the local curvature, the 
reciprocal of the radius of the circle passing through the wake point in 
question and the two wake points adjacent to it. 

The basic equation for velocity computation of the main wako (tip 
vortices) , inner wake, bound vortices and image wake, at an arbitrary point 
P was taken as 

V= r (Cos e -Cos e )/(4 h) 

1 2 


and for self-induced velocity to be 


V= r (Ln{(8H'a )*Tan{ 4> /4)-l.J/(8 n R) + 
111 

y-’[Ln((SV'a )*Tan( ^ /4)-l.]/(8:r r) 
2 2 2 


Ihese two equations will be fonnulated in three dimensional space in terms 
of coordinates of the segnvents and their vorticities, so the velocity com- 
ponents induced hy different segments at a point can be summed, 

Li this chapter all necessary- formulas required by the computer code 
will be derived and all the flight parameters will be calculated. 


3-1 INDUCED VELOCITY OF A VORTEX SEGMENT AT POINT P IN OHREE DIMENSIONS. 

In equation (2-4) the velocity vector induced by segment S at the ’ 
point P is perpendicular to a plane containing the segment and point P. As 
we go from one segment to another the direction of the plane changes, conse- 
quently the velocity vector will change direction too. In order to 
cal iulate the total velocity at point p, it is necessary to formulate the 
components of velocity vector. The components of the total velocity is the 
sum of the components of induced velocities of all the segments at p. 

In three dimensional space, point P is defined as P(x,y,z). Ihe two 

ends of segment S, are points P and P with coordinates (x ,y ,z ) and 

12 111 

(X ,y ,z ) respectively. 

2 2 2 


The velocity vector V at point p is perpendicular to the plane of P, 

P and P . Ihe unit vector in velocity direction may be calculated as: 

1 2 


P =(r XL)/ 
1 


r X L 
1 


(3-1) 


also 

|r X l| = r L sin 8 = h r I (3-2) 

' 1 ' 1 1 1 ■ 

assuming 
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A ^ 

TxlT = w i+ 

1 12 3 

Ihe unit vector can be decomposed as 

P = y /(L h) k=l,2,3 (3-3) 

k k 

V , V and 1 / are the components of vector (Tx L ) in X, Y and z 
12 3 1 

directions. Using equation (3-3) in (2-4 ) , components of induced-velocity 

of a segment at point p may be written as 

2 

V = r (Cos e -K:os e ) l' /(4 7T L h ) k=l,2,3 (3-4) 

k 1 2 k ^ 

where Cos 0 +Cos $ and h may be C£dculated in terms of r , r and L! 
12 12 

2 2 

Cos e 4Cosd = (r +r ) [L -(r -r ) ]/(2 r r L) (3-5) 
121212 12 

2 2 2 2 2 2 

h = [(r +r ) -L ] [L -(r -r ) ]/(4L ) (3-€) 

12 12 

By substituting (3-5) and (3-6) in (3-4) components of induced velocity can 

be computed in terms of coordinates of points p, p , p and t. 

1 2 

2 2 

V = r (r +r )/[27rr r ((r +r ) -L )] y k=l,2,3 (3-7) 

k 12 1212 k 

where 

2 2 2 2 
r =(x-x ) +(y-y ) +(z-z ) 

1111 

2 2 2 2 
r =(x-x ) +(y-y ) +(z-z ) 

2 2 2 2 

2 2 2 2 
L =(x -X ) +(y -y ) +(z -z ) 

12 12 12 
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and 




F'. 


, 4*1 


OF 




V =(y-y ) (z -z )-(z-z ) (y --y ) 

1 112 112 

V «(z-z ) (X -X )-(x-x ) (z -z ) 

2 112 112 

V »(x-x ) (y -y )-(y-y ) (x -x ) 

3 112 112 

3-2 SELF-INDUCED VSLOCm IN raREE DI!®ISION 

In equation (2-6) the velocity vector V is perpendicular to the plane 

containing segment S and S . The unit vector in velocity direction can be 

1 2 

computed similar to the method given in the previous section. 


A 

P =(rx L )/lL X L I (3-8) 

12 12 

If 

^ =(y -y ) (Z -Z )-(z -z )'{y -y ) 

1 12 23 12 23 

V =(z -z ) (x -X )-(x -X ) (z -z ) 

2 12 23 12 23 

P =(x -X ) (y -y )-(y -y ) (x -x ) 

3 12 23 12 23 

then 


- - A A A 

L X L = p j+ i/ k' 

12 12 3 

and 

22 2 1/2 

p = /[V +V +P 1 (3-9) 

k k 1 2 3 

the radius of the circle formed ty s and S may be calculated from 

1 2 

following formula. 
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2 2 2 2 

R =(L *L *L )/[ (L + L - L ) (L +L +L ) (-L +L +L ) 

123 123123 123 

(L -L +L )] (3-10) 

12 3 

where 

2 2 2 2 

L =(x -X ) +(y -y ) +(z -z ) 

112 12 1 1> 

2 2 2 2 

L =(x -X ) +(y -y ) +(z -z ) 

2 2 3 2 3 2 3 

2 2 2 2 

L =(x -X ) +(y -y ) +(z -z ) 

3 13 13 13 

<p and d> can be computed from 
1 2 

-1 

=2 Sin [L /(2R)] (3-11) 

1 1 

-1 

0=2 Sin [L/(2R)J (3-12) 

2 2 

Substitution of equations (3-10) to (3-12) in (2-6) provides the com- 
ponents . of the self-induced velocity of an arc oontaing S and S in terms 

1 2 

ot r t r f coordinates of p , p and p . 

12 12 3 

3-3 VftRIATION OF dROJLATION Wim AZIMjm ANGLE. 

In hover with no wind, the velocity of the air with respect to a point 

on a blade remains constant with azimuth angle: 

V{xlr,t)=Qz 

Assuming that the induced velocity at the rotor disc is snail in comparison 
with the velocity of the blade, Kutta-Joukowsky law may be used to calculate 
elementary thrust and elementary rolling manent; 


dT= P r V dr 
dM »P r" V r Sin^dr 

X 

Integration of these two equations over the entire disc (with the assumption 

of constant vorticity along the blades and the fact that the thrust ap>- 

poximately equals the weight of the aircraft) yeilds; 

r=(2W)/(bPV R) (3-13) 

.t.. 

and 

M =0 

X average 

Here b is the number of the blades, 

V is the velocity of the tip of the blades and 
t 

R is the rotor radius. 

In the absence of the wind or in hover equation (3-13) may be used for 
computation of circulation F . But in the presence of the wind or in for- 
ward flight it may not be used, because the velocity of the air relative to 
a point on the blades is not constant and varies with azimuth angle. 

V(i/;,r)= r+v iJSinip (3-14) 

f 

Here V is the forward speed or wind velocity. Ihere will be some average 
f 

rolling moment per revolution if it is assumed that the vorticity does not 
vary with azimuth angl<s. A^oiroetry in velocity profile with constant p is 
the cause of the appearance of the rolling moment. Ihe rolling moment can 
be found ty integrating the infinitesimal mcment over the entire disc. 

dT(0,r) =pr (i2r+V sinxH) dr 

f 

dM ,t)=Pr (i2r+v Sirup)sin \}j rdr 

X 
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Integration along the blade at a constant azimuth angle provides; 

2 

T(\i») =(Pr R )(l+2uSin\if)/2 (3-15) 

3 

M (tp)«{P r R )(l+1.5A<Sin^^)Sini/; /3 (3-16) 

X 


where 

» V /( R) 
f 

Finally the average rolling monent per revoltion 

2 

M =-bPrRV/4 (3-17) 

X average t 

Thrust offset can be found by eliminating r between the thrust and the 
average rolling moment. 

y =0.5 P 
0 

As forward speed or velocity of the wind increases the thrust offset 
increases too. Such large thrust offsets are not realistic for conventional 
helicopters, although they might be tolerated in seme other configurations 
like side-l3y-side helicopters. 

Elimination of average rolling moment can be achieved by postulating 
that the blade thrust mcment with respect to its flapping axis remains 
constant. 

3 

M =(Pi? R /3) r(l+1.5psin ip- )=const. 
f 


Ihis requires that circulation to be 

r=(const.)/(l+1.5PSin^P) (3-18) 

Equation (3-18) have been employed in the computer code for blade cir 


eolation at different ijr . 


\ 


3-4 CORE SIZE 

As a point approaches the center of a potential vortexr the velocity 
induced at that point fcy the vortex tends to go to infinity. Because the 
velocity of the air is finite and can not reach infinity, this tendenc/ is 
not realistic. Itierefore some modifications are needed to overcome this 
problem. In chapter 2 the Scully model [3] 

2 2 
r =r(Va) / u. + (Va) ] 
c 

was assumed to have proper characteristics. This model required the 
knowledge of core radius and intensity of the vortices. 

In reference [2] an average value for radius of fully rolled up tip 
vortices has been suggested as 

a=0.003 R (3-19) 

The results were obtained by experimental data and verified to be 
quite accurate for high aspect ratio rotor blades, in the present stu<^, 
rotors have been assumed to have high aspect ratio blades, the tip vortices 
are assumed fully rolled up from the time they are generated, and the ground 
has no effect on the core size. With these assumptions the average value 
of a=.003 R was used as starting value for the vortex segments when they are 
at the tip of the blades. The core radii of the vortex segments at other 
places are calculated by assuming that the volume of a vortex segment 
ranains constant as time passes. Therefore 
2 2 

a = a L / L (3-20) 

2 112 

where L is the length of the segments. 

Considering the fact that there is no air where the blades are and the 
bound vortices are blade replacements, the core radius of the bound vortices 
was chosen to be one-half of the chord length. 
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unlike the tip vortices or bound vortices, the inner wake is not a 
single concentrated vortex line but a continous sheet of vortices which was 
replaced by a vortex filament, Hie velocity induced by this filament at 
neighboring points is quite high considering the smooth velocity profile in- 
side the wake. To overcome the unrealistic behavior around the center of 
this filament, the radius for core of these segments was chosen to be 
a«0,03R. A few other core sizes for. inner wake were examined. The results 
have indicated the choice of the analler radius will result in numerical in- 
stability which is one of the main difficulties of this approach. 

3-5 INCREASE IN ACCURACY AND PASTER ALGORITHM 

Another major difficulty of free wake approach is the consumption of 
a large amount of computation time. If N is the total number of the seg- 
ments constructing the wake, the computation of velocity at only one point 
including the image wake requires 2N times ciilculation of equations (2-4) or 
(2-6) which themselves require number of operations. There are N points on 
the wake whose velocities are computed at each time slice. Therefore at 

2 

each iteration 2N times equations (2-4) or (2-6) are computed. 

Considering time for integration of velocities and other procedures, com- 

3 

putation time becomes proportional to N . 

On the other hand, the accuracy of the calculations greatly depend on 
the number of the segments per revolution and number of the revolution in 
the wake. The shorter the segment length the better the results. Therefore 
there is a trade-off between accuracy and computation time. 

Most of the cxjntribution of induced velocity at a point comes from the 
segments in the vicinity of that point. To increase the accuracy of these 
velocity contributions, each segment close to the point of interest has been 
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broken into anedler segments and velocity induced ty each smaller segnent 
has been completed and added together. 

The segment division was done fcy passing a circle through two neigh- 
boring segments and dividing the arc passing through the end of each segment 
into two equal parts. 

Also if h/R >2.5, the contribution of induced velocity of this seg- 
ment in comparison with the contribution of another segment with VR < 0 \ 
is negligible. Therefore it is not necessary to compute the velocity con- 
tributions of segments with distances beyound ^1^2, 5. 

3-6 COORDINATES OF THE IMflGE ?7AKE 

As previously mentioned, to include the effect of the ground on 
helicopter rotor aeroctynamics a mirror image for the wake has to be assumed. 
The velocity induced by this wake is computed the same way as for the main 
wake with the exception of the self-induced velocity. The image wake was 
broken into . segments exactly like the main wake, so that segments of the 
image wake are the mirror image of the segments of the main wake. 

Assuming the rotor can only tilt forward, the coordinates of the mir- 
ror image of point P(x,y,z) for the rotor with (Tip Path Plane angle) and 
H (the height above the ground) can be obtained from the following 
equations. 

X = X Cos(2a ) - z Sin(2a ) - 2 H sin or (3-21) 
mi t t t 

y = y (3-22) 

mi 

z = - X Sin(2a ) - z Cos(2 0i ) -- 2 H Cos a (3-23) 
mi t t t 

X , y and z are the coordinates of the mirror image of point p. 
mi mi mi 

3-7 r ASSIGNMENT 


- 35 - 




ORIQWAt PAGE IS 
OF POOR QUAUTV 


Theoretical studies of helicopter blades have indicated that there is 
a direct correlation between load of the rotor and the circulation of the 
blades. If for simplisity circulation is assumed to be constant along the 
blade, equations (3-13) or (3-15) can be used to calculate average 
circulation. For a two-blade rotor 


2 

r /(i?R ) C (3-24) 

m T , 

where cis defined from 

T 

2 2 

C = W /[(7TR) P (7TR )] 

T 


Scmewhat better results are obtained if elliptic profile is assumed 
for spanwise blade circulation. 


2 

r /( R ) = 4 C (3-25) 

m T 

However even more realistic and better results can be obtained using 

experimental data. Figure 3-3 was plotted using data extracted from 

reference [2] . The slope of a nearest straight line to the points extracted 

from experiments is assumed to be a better relation between C and p . 

T m 


2 

r /(i?R ) = 5.075 C 
' m T 

Considering the variation with asimuth angle, we obtain 


(3-26) 


2 

r /(S3R ) = 5.075 C /a + 1,5 P sin^^r) (3-27) 
m T 

Equation (3-26) has been implemented to assign circulation to the vortices 

released at the tip of the blade at the time of generation. 

To be able to compute the vorticity of each segment at a later time, 

the vorticity dissipation should be known. Since the total time which a 

vortex segment is in the domain of velocity computation is less than two 

seconds, the dissipation of vorticity can be neglected. The assumption of 
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constant vorticity provides a simple formula for vorticity computation of 

the segments as they shrink or elongate. 

p =r L / L (3-28) 

1 2 2 1 

3-8 CXX)RDINATE DESOilPTION AND NaffiNCLATORE 

A coordinate system fixed in the tip path plane of the rotor is used. 

The model for a two blade rotor and its wake is shewn in figure 3-4. As 

noted on the figure a free stream of dimensionless magnitude fi is directed 

at an angle a to the tip path plane. The azimuth angle position ip of 
t 

the rotor is defined to be the angle between blade vortex number one and X- 
axis, as shown in the figure. The point P is the wake reference point; 

ij 

the first subscript, i, increases successively proceeding down the wake vor- 
tex for a given blade, and the second subscript, j, denotes the blade number 
of the blade which generates that wake vortex. Each wake segment is 
sociated with that end point having lower first subscript. In the 
coordinate systan described above the location of point P is then defined 

ij 

ty X / y and Z . The total velocities associated with the X, Y and Z 
ij ij ij 

directions are defined U , V and V? respectively. And finally each 

ij ij ij 

element (i,j) is assigned a dimensionless core radius a and strength r . 

ij ij 

3-9 NUMERICAL DAMPING 

The goal in the present approach is to compute the location of the tip 
vortices ty iteratively computing the velocities and the new locations of 
the tip vortices. The computation is terminated when a periodically steady 
state solution for the locations of tip vortices is achieved. 

Succesive computation of velocities and locations does not always 
result in a steac^ state solution. If flight conditions are such that in 


- 3 ^ 



reality the wake is unstea<^ or if interaction of the old vortices and newly 
generated are present, then the computations will result in an unstable 
solution which does not necessarily represent the true answer. It is also 
possible that the wake is stable in reality, but its numerical iteration 
is unstable due to either method of integration, or poorly chosen initial 
conditions. 

To overcome these problems, .and also to be able to obtain an average 
location for the case in which the wake is unstable in reality, numerical 
damping was introduced. 

Detailed study of the main wake indicated that the computation in- 
stability was started and magnified when a vortex segment of the wake of the 
blade number one and a vortex segment of the wake of the blade number two 
moved very close to each other. Because the velocity induced by a vortex is 
inversely proporti.onal to the distance from the vortex, the two vortices in- 
duce large velocities on each other. In the next step of integration they 
unrealistically move far from each other. This behavior may deform the wake 
such that the continuation of the computation only worsens the results. It 
is possible to introduce an upper limit for the velocity to prevent large 
wake deformation. This can stabilize the iteration but convergence becomes 
very slew, therefore this concept was rejected. 

In steady state there should be no difference between each successive 
blade. If at a particular time when blade number one is at some azimuth 

angle, point p on the wake blade number one is at point (x , y , z ) , then 
Ij 0 0 0 

when blade number two is at the same azimuth angle, point p should be at 

2j 

the same point (x , y , z ) . Based on this criteria numerical damping can 

0 0 0 

be introduced. In attanpts to reach faster steady state solutions, if the 


location of a point on wake of blade number one was too far from the loca- 
tion of the corresponding point for the wake of blade number two, a point 
between these two locations was used as a better estimate of the new loca- 
tion of the two points. 

This method was employed in our computer program. Use of a point 
from 65% to 75% of the distance between the old location and the newlocat- 
ion gave the best convergence and mo^ stable solution. 

3-10 SIKJCrURE OP THE OOMH3TER PROGRAM 

The computer code documented in appendix B has been constructed to 
simulate the physical flow fcy means of the modeling the main rotor, the tip 
vortices and the ground. Given an initial wake geometry and aircraft flight 
condition, it proceeds to compute the velocity and integrate in time suc- 
cesively until a periodically steady solution is obtained. 

The program has been written so that it can be used for out-of -ground- 
effect (CXsE) , as well as in-ground-effect (IGE) . By proper usage of flags, 
the wake can be shortened or enlarged; the iteration can be applied to a 
portion of the wake or to the whole wake. The flow of information as com- 
putation proceeds, is presented schematically in figure 3^5. 

At the beginning the program reads the flags and flight conditions. 
If an initial wake is available then the program reads the wake and its 
properties; otherwise it creates a simple helix for use as an initial wake. 
Subroutine "VELOC" computes the total velocity at the vortex segment's ends; 
subroutine "NEl'JLOC" computes the new location of the end points and sub- 
routine "SMDOTH" (which is called fcy "ND-njOC" subroutine) applies the 
numrical damping. If an acceptable solution is reached the program is ter- 
minated, otherwise it continues until the number of iterations exceeds a 
given maximum. 
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THE MAIN PROGRAM 



FIGURE 3-5. Flow-Chart Of *nie Conputer Program 



The description of what the subroutines do as well as their inputs and 
their outputs are given at the beginning of each subroutine as comment 
statements. The main program and all subroutines are documented in appendix 


B. 

3-11 FLIGHT OONDITIDN PARAMETERS 

The formulation of the model was nondimensionalized for the purpose of 

coding, with lengths made dimensionless ty rotor radius "R" and velocities 

ty rotor tip speed "V The flight conditions of the aircraft being repre- 

t 

sented relate to computer program through the following paraneters: 

/^ = V /V advance ratio 

f t 

2 


r = r / (i?R ) 

vorticity 

a 

tip-path-plane angle 

t 

a/r 

core radii 

2 4 

C = /( PP r R ) 

thrust coefficient 

T 

H/R 

height above the ground i 

R/C 

aspect ratio. 
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OiAPTER 4 


AERQDYMAMIC RESULTS 

Hie computing program given in appendix B, computes the location, 
velocity, core radius and circulation of the fully rolled up tip vortices as 
they are generated and moved in space. Later modifications to the progrem 
all wed the computation of the velocity of any particle in space, and as a 
consequence its new location after a short period of time. This enabled the 
plotting of streamlines or streaklines. Some of the results of this chapter 
were obtained using the modified version of the progran. 

It should be remoiibered that this program was originally developed for 
stability derivative computation in ground effect and not for detailed stucfy 
of the aeroc^namics of the rotor. Therefore, a number of assumptions oon^ 
sistent cur purpose was made. The assumptions were: 

a. potential flew, 

b. effect of fuselage and tail rotor were neglected, 

c. tip vortices are fully rolled up as they are generated, 

d. circulation along the blades changes in two steps, 

1. at the tip r / R = 1.0 and 

2. at the location r / R = 0.7 

e. blades were replaced by line vortices, 

f. no decay in circulation and no merging of vortices, 

g. vortex segments are assumed to be straight lines and 

h. the wake is stable. 

4-1 UNSTABLE WAKE (^D NUMERICAL DAMPING) 

As was mentioned in the previous chapter, with no numerical damping an 
unstable wake might result from successive computation of the velocities and 
locations, in figure 4-1 the cross-section of the location of the main wake 
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— — Iteration Cycle 1 

-~A-— — Iteration Cycle 2 
• •0 Iteration cycle 3 


Figure 4-1. Numerically Unstable Wake For Cycle 1, 2 and 3 
(No Numerical Damping) 




in hover cuid out of ground effect for cycles 1, 2, and 3 of iteration with 
no numerical damping are shown. As can be seen, in iteration number 2 the 
end points of two segments of the wakes of the two blades have approached 
each other closely. Because point (2,13) is close to segment (1,7), the to- 
tal induced velocity at this point has wrong magnitude and is in the wrong 
direction. In the next cycle it circles the segment (1,7) (see figrue 4-1) . 
As a consequence the whole wake is distorted and if the computations con- 
tinue similar instabilities are repeated. In ground effect the results were 
worse because of additional interaction of the old and new vortices. 

4-2 STABLE WAKE 

The use of numerical damping resulted in a stable wake for a rotor in 
the same cjonditions as the previous section. Figure 4-2 and 4-3 show the 
wake cross-section and three-dimensional wake geonetry after reaching a pe- 
riodic stea<^ state solution. Figure 4-4 is the same rotor in ground effect 
and in hover, for this case again a periodic steac^ state solution was 
obtained. 

4-3 COMPARISON OP THEORY AND MEASUREMENT 

To check the results obtained from the computer code, they were com- 
pared to the results obtained from wind tunnel tests performed in reference 
[ 1 ]. 

As can be seen in figure 4-5 the computational results are very close 
to the measured ones, especially in the first cycle which is the area of 
interest. The error between the new wake and the measured wake increases as 
one goes downstrecm. The reason for this discrepancy^ is the error in for- 
ward integration and number of the assumptions made. However this error is 
not very important, because the effect of the far vortices on the rotor disc 
is small compared to that of the newest cycle. 
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Figure 4-3. Three Dimentional wake Out Of Ground Effect 
(0^=0.0037, a^=0.0, Numerical Damping=70%, 
And Aspect Ratio=13.7) 










I^ieoretical 




Measured 



since experimental c3ata for the main wake in ground effect were not 
available, quantitative comparison in ground effect was not possible. But 
the follcwing results indicate there is a good comparison of the theory and 
the actual flow patterns. 

4-4 PLa; FIELD IN GROUND EFFECT PND HWER 

Experiments with a radio controlled helicopter shewed the snoke coming 
from the engine exhaust would follow different patterns in different flight 
conditions. 

For hover and very close to the ground, the smoke flowed upward near 
the center of the the rotor and down ward elsewhere inside the disc area. 
As the rotor moved further away from the ground, the amount of smoke going 
upward around the center was reduced. For h / R > 1.0 no smoke flowed out 
of the center of the rotor, as can be seen in figures 4-6 and 4-7. 
Theoretical results show the same behavior, as can be seen in figures 4-8 and 
4-S. 

4-5 GROJND VORTEX 

Experimental studies of a rotor in ground effect and in hover have 
shown that the vortices generated by the tip of the blades close to the 
ground wil] move away from the rotcr, and there is no vortex interaction and 
no air circulation through the rotor. But in the presence of wind or in 
slow forward flight, the vortices close to the ground upstream of the rotor 
will roll up and form a horseshoe vortex around the rotor. This is called a 
ground vortex and can create problans if there is gusty wind. 

For slew forward speed or light winds, the ground vortex will stay far 
upstream. As the wind velocity increases, it moves closer to the rotor. 
Interaction between old vortices and new vortices increases as forward 
speed increases. At a certain speed (depending on the height above the 


ORlQf^’ril PAC>' '9 
W POOR QOAUVf 



Figure 4-6. FICV Visualization m 


Ground Effect i;/R-o.4 



Figure 4-7. Flow Visualization in Ground Effect H/R>l.o 
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ground and tip-path-plane angle ) , when the ground vortex is nearly under the 
leading edge of the rotor, the interaction reaches its maxiituan. If the for- 
ward speed increases further, the ground vortex will be washed aw^ down- 
stream. Figures 4-10 and 4-11 show the actual flew at different speeds 
(reference [5]). Similar results were obtained by the present approach. 
These results are shown in figures 4-12 and 4-13. 

4-6 CONEPCmCE 

Detailed stu(ty of the locations of the points on the main wake with no 
numerical damping showed two interesting results which resulted in a method 
of applying numerical damping which was described in the previous chapter. 

First, the points closer to the rotor disc converge to their steac^r 
value faster than the point far from the disc. Secondly, the convergence to 
the final values were oscillatory. Figure 4-14 shows the convergence of the 
wake after applying numerical damping. As can be seen, the points closer to 
the rotor converge faster to their final values. It was also concluded that 
recomputation of the points which have reached their final location is a 
waste of time. Therefore, another check was added to the progran to skip 
the points whose position changes were very small. This made the program to 
run faster for each flight condition. 
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Figure 4-10. Ground Vortex Visualization (FyR»0.53, AR=9.8 And 
A =0.05 Pefcrence [5]) 



Figure 4-11. Ground Vortex Visualization (H/R=0.53, AR=9.8 And 
A =0.055 Reference [5]) 
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Figure 4-14. Wake Convergence Out Of Ground Effect With 70% 
Numerical Damping 



CHAPTER 5 


onNCT.usinyts 

Ihe model and the computer program developed in this stu^ provides 
the velocity, location, and circulation of the tip vortices of a two-blade 
helicopter in and out of the ground effect. Comparison of the theoretical 
results with seme experimental measurenents for the location of the wake in- 
dicate that there is excellent accuracy in the vicinity of the rotor and 
fair amount of accuracy far from it. Having the location of the wake at all 
times enables us to compute the history of the velocity and the location of 
any point in the flow. The main goal of our stucty, induced velocity at the 
rotor, can also be calculated in addition to stream lines and streak lines. 
Since the wake location close to the rotor is known more accurately than at 
other places, the calculated induced velocity over the disc should be a good 
estimate of the real induced velocity, with the exception of the blade loca- 
tion, because each blade was replaced only by a vortex line. 

Because no experimental m.easuranents of the wake close to the ground 
were available to us, quantitative evaluation of the theoretical wake was 
not possible. But qualitatively we have been able to show excellent 
agreement. Comparison of flow visualization with our results has indicated 
the location of the ground vortex is estimated excellently. Also, the flew 
field in hover is well represented. The addition of numerical damping 
provided the three important following results; 

1. faster convergence, 

2. stable numerical solutions for stea^ wake and 

3. computation a tim?. average location for an unstable wake. 

These results for stable and unstable wakes in steaefy flight conditions 
should be accurate enough for stability derivative computations. 
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The computation of stability derivatives requires the knowledge of the 
aeroc^namic forces and moments around a trim condition. As a consequence 
induced velocity around the trim condition should be known. Therefore# the 
computation of stability derivatives requires the execution of the present 
computer code for a variety of different flight parameters around the trim 
condition. As a result# that portion of stu(^ requires a considerable 
amount of computer time# the lack of which has delayed progress on computa- 
tion of stability derivatives. 

As mentioned in previous chapters# because of the importance of the 
old vortices in low altitude and low advance ratio# and also because of the 
consideration of the mirror image wake the number of points in the wake 
should be enlarged. The time required for each iteration was proportional 
to the cube of the number of points in the wake. Although the progrem has 
been made to have faster convergence# it still requires a large amount of 
computing time for near hover conditions. However for high advance ratios# 
or for high altitude the vortices do not interact with each other or with 
the rotor#* therefore a smaller wake is sufficient for stability 
computations. 

This program can also be employed for two-blade propeller stuc^ in or 
out of the ground effect. With some modification# it can also be used for 
rotors and propellers with more than two blades. A fair amount of modifica- 
tions is required to employ this progran for the interaction of two rotors 
or two propellers or more. The modifications are only in usage of different 
subprograms and m.emories required, not the formulation or method of 
approach. 
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The preliminary stu<^ of a cylinder vortex sheet in ground effect has 
indicated that when a rotor in hover is close to the ground, the ground be- 
haves like a spring with a damper. The spring constant increases, as the 
rotor approaches the ground. A better understanding of the ground effect in 
hover and forward flight can be made with the help of the present program. 

A quantitative validation of the present stu^ could have been 
carried out if more experimental data of measured wake in ground effect were 
available. However for qualitative stu<^, the present work is a very good 
tool for prediction of the ground vortex, and computation of the induced 
velocity. 
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APPENDIX A 


OF POOR 


A-1 SELF-INDUCED VaOCI'W 

The self-induced velocity of segment APC in figure a-1 may be computed 
by subtraction of the induced velocity of segment ABC at point P from the 
self-induced velocity of a whole vortex ring with the same curature as seg- 
ment APC. 

The Biot-Savart Law for the induced velocity of segment ABC at point p 
can be caculated as follows: 

^0 _ 3 

(d X S) / d (a-1) 



Figure A-1. Vortex Ring 
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OF POOR 8 '' 


Assuming point P and the whole ring is in X-Y plane, then 
^ = R (Cos^-1.) 1 + R Sin<^ j (a-2) 

ds = R t - Sin0 i + Cos^ jj (a-3) 

And 


3 3 3/2 

d e- R [2-2 CO30 ] 


(a-4) 


Using equations (a-2) through (a-4).. in the induced velocity equation (a-1) 
can be reformulated as an elliptic integral. 

1/2 


V = r /(2 
P 

Let 20 then 


Jbttr) 

\ 


d«^ /[I. - Cos^J 


2 1/2 


V = r/(4 7T d0/[l. - Cos<^ ] 

i/2 


(a-5) 


(a-€) 


or 


V = - r /(4 R) Ln( Tan( /4)) 
P 


(a-7) 


Subtraction of V from self-induced velocity of the whole vortex ring 
P 

V =r/(4n-R) [Ln(8R/a)-l.) (a-8) 

s 

leads to calculation of the self-induced velocity of segments APC. 

V = r/(4 7r R)[Ln(8EV'a Tan( 0^/4))-!.] (a-9) 

P 
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APPENDIX B 


B-1 OPERATIONAL INFORMATION EOR TOE COMHJTER PROGRAM 

program presented in this appendix was written in Fortran IV and 
executed on DEC-20, IBM-370, and VAX-11. 

The inputs to the program are as follows: 

^XX)MP Number of point to be iterated on 
IC Number of blades 

Hi Number of segments for each blade wake 

NOONS Number of points that does not need tb,e computations 

lADD Set to 0 if wake should not be enlarged 

MIRIM3 Set to 0 for in ground effect ans 1 for out of ground effect 

PSIO Initial azimuth angla 

REV Number of revolution to be iterated 

CT Thrust Coefficient 

XMU Advance ratio 

ALFHAT Tip-path-plane angle 

RB Aspect ratio 

H Height ABove the Ground 

ALFR Damping percentage in r direction 

ALF2 Damping percentage in z direction 
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o-0ct-l?83 02:S3:09 %Mi FORTRAN V3,i-23 
2£-Jul-1933 23:02!4? SYSIUSER!CEABERI,F,EP3KAIH.F0RI8 


0001 

0002 

0003 

0004 

0005 
OOOi 
0007 
0003 

0009 

0010 
0011 
0012 

0013 

0014 

0015 
OOlo 
0017 
0013 

0019 

0020 
0021 
0022 

0023 

0024 

0025 
C02i 
0027 
0029 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 


C THIS IS THE LAST VEEION OF KELLICAL WAKE TILL JUNE 83, 

C F0R39 IS INPUT FILE READING FILE (IKPUT.DAT), 

C F0R40 IS STORING DATA FOR NEXT STEP (INTER, DAT), 

C F0R50 IS FOR BLADE VORTICES LOCATOINS (OUTPUT), 

c mm HAiN PROGRAN um 

C 

C THIS PROGRAH CALCULATES THE INDUCED VELOCITY OF A HELICOPTER 
C ROTOR IN GROUND EFFECT, 

COHHON/ALLSUB/ALFO»ALPHAT»CAT»CAT2»CA2TfCTjDALFO»DF’SI»EPSr 

1 H*IADDfIAV>HIRIKGiNA)NCONSrNDPSI*NPP.>NRyfNU>NUhPh 

2 PSIfP3IF.PSI0)RAD,RC»REV)SAT,3AT2»SA2T»TKCATr 

3 TFRCiTHSATFXMU»XX»YYFZZfALFR»BETR»ALFZfBETZ» 

4 A(180;2)FGAHA(180F2)FGAHB(30)FSEG(180r2')FU(180F2)f 

5 y(lS0F2)FW{180F2)FX(lS0F2)FY(180F2)FZ(I80F2) 
COHKON/MIS/ SX(180f18)fSY(130f13)fSZ(180f18) 
COHHOH/NTURN/NCOflP 

PI=4,«ATAN(1.) 

RAD=PI/180, 

NCT=0 

C 

C READS ALL THE FLAGS fFLIGHT CONDITIONS fNUNERICAL 

C DAHPING PROPERTIES fCOR OF EACH SEGMENT OF THE 

C CIRCULATION OF EACH SEGMENT AND HAKE LOCATION 

C FROM PREVIOUS RUN. 

C 

READ(39f») HCOMP 

READ ( 39 F I) NA F NH F NCONS f IADD fMIRIMG fNDPSI f NPR 
READ(29fI>PSI0fREVfCTfXMUfALPHATfRBfH 
READ(39f*)ALFRfBETRfALFZfBETZfEPS 
URITE(50f1000)NAfNWfNC0NSfIADDfMIRIMGfNDPSIfNPP. 

HRITE(50f1004) FSIOfREVfCTfXMUfALPHATfRBfH 

URITE(50Fl00i)ALFRFBETRFALFZFBETZFEPS 

HW1=NH41 

NRU=NH/NA 

RC=2.*RB 

101 IF(PSI0,LT,1) GO TO 301 


0033 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 
0017 
0043 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


READ(39Fl)((A(lFJ)FiaFNHl)FJ=lF2) 

READ(39f»)((GAMA(IfJ)fI=1fNH)fJ=1f2) 

READ(39fI)((X(IfJ)fY(IfJ)fZ(IfJ)fI=1fNH1)fJ=1f2) 

READ(39Fl)((SX(lFJ)FSY(lFJ)FSZ(lFJ)Fl=lFNUi)Fj=lFNA) 

NAHP=NA/241 

IF(IADD.EQ.O) GO TO 301 
CALL ADDP(NAHP) 

HH=NH1-1 

A(NU1f1)=A(NHf1) 

A(NH1f2)=A(NHf2) 

GAMA(NWf1)=GAMA(NH-1f1) 

GAMA(NWf2)=GAMA(NW-1f2) 

301 CONTINUE 
[iPSI=2,»PI/NA 


—■ VORTICITY OF THE FIRST SEGMENT VERSUS AZIMUTH OF THE ROTOR — - 
BLADE. 


OF POOK QU'i«-riY 


SIE=0. 

DO 299 1=1 fNA 
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OF POOR QUALITY, 


6-0ct-1933 02!53',09 WX-ll FORTRAN V3.1-23 
26-JUM9S3 23!02.N9 SYSfUSERIOSABERI.REPJHAIN.FORJS 


0058 

0059 

0040 

0041 

0042 

0043 

0044 

0045 
0044 

0047 

0048 

0049 

0070 

0071 

0072 

0073 

0074 

0075 
0074 
0077 
0073 
0079 
0030 
0081 

0032 
0083 

0034 

0035 
0034 
0037 

0033 
0039 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 
0093 

0099 

0100 
0101 
0102 

0103 

0104 

0105 
0104 

0107 

0108 
010 ? 
0110 
0111 
0112 

0113 

0114 


S1=SIN(SIE) 

S2=SIH(SIE4DPSI) 

GAH=,5/( 1 .41 .5»XNU»Sm. 5/(1 .+1 . 5*XHU*S2) 
GAHB(I)=5.1IGAH/(2,»P!) 

SIE=SIE+DPSI 
299 CONTINUE 

SAT=SIN(ALPHAT»RAD) 

CAT=C0S(ALPHATIRAH) 

SAT2=SAT*SAT 

CAT2«CAT»CAT 

SA2T=SIN(2.tALPHAT*RAD) 

CA2T=C0S(2.«LFKAT*PAD) 

THSAT=2.mSAT 

THCAT=2.»HICAT 

TMP=RC»IIPSI 

THPl=SDRT(THP»(THP+2.)) 

TFRC={ THP-TKPl +AL0G ( 1 . 4TNP4TNP 1 ) ) /DPSI 
PSIFF=PSI0+340.tP,E'.' 

PSI=PSI0tRAH 

PSI0=P3I 

THE FIR3T TIME OF A FLIGHT CONIUTICN INITIAL WAKE IS 

CALLED TO PRODUCE AN INITIAL CONDITION AN INITIAL 

CONDITION FOR THE WAKE. 


PSIF=PSIFF»RAD+.05 
15 IF(PSIO.LT.O.Ol) CALL INWAKE 

COMPUTES LENGHTH OF EACH SEGMENT. 

DO 250 J=l>2 
DO 250 I=l»NW 
11 = 1+1 

DX=X(IhJ)-X(I>J) 

DY=y(IljJ)-Y(IfJ) 

DZ=Z(I1»J)-Z{I»J) 

SEG ( I > J ) =SQRT ( DXIDXtDYID Y+DZ tDZ ) 

250 CONTINUE 

START OF THE MAIN LOOP OF THE MAIN PROGRAM. 

SUBROUTINE VELOCITY COMPUTES THE TOTAL INDUCED 

VELOCITY AT ALL THE POINTS OF THE HAKE 

20 CALL VELOC 
40 IF(NCT.NE.O) 60 TO 50 

TN=IFIX((PSIt.05)/DPSI)/(1.»NA) 

WRITE(50fl002)TN 

NAH=NA/2 

WRITES THE CROSS SECTION OF THE WAKE 

PLANE X-Z 

NWE=NHHIAH 
DO 10 I=lfNHE»HA 

URITE(30.1005)l7X(I»l)»Y(Irl)»Z(Ifl)fU(Irl) 

K=I+NAH 

10 URITE(5v;1005)K»X(Kf2);Y(Kf2)»Z(K»2)FW(K»2) 


~6<h- 
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6-0cl-1933 C'2|SS:09 VAX-11 FORTRAN V3.1-23 
26-Jul-19S3 23;02*.-l9 SYSJUSERICSAEERI.REPTHAIN.FORIS 


0115 «RITE(50»1005)NMlfXaiMl>l)»Y(NUl»l)jZ(Hyirl) 

0116 «RITE(S0j1000) 

0117 DO 30 I=l»NOEfNA 

0118 aRITE{50>1005)I»X(I;2)>Y(I»2)iZ(I»2)»U(I»2) 

0119 K=HNAH 

0120 30 URITE'.50»1005)KiX(KFl)»Y(K»l)fZ(K.l)f0(K»l) 

0121 URITE(50rl00S)NUlFX(NUlr2)iY(NUlr2)FZ(NaiF2) 

0122 URIT£(50»1000) 

0123 50 NCT=NCT+1 

012< IF(NCT.GE.NDPSI)NCT=0 

0125 PSI=PSHDPSI 

0126 IF(PSI.LT.PSIF)GO TO 30 

0127 C 

0128 C STORE PROPERTIES OF THE fc'AKE ON A FILE TB-SE 

0129 c TO BE READ FOP. THE NEXT RUN, 

0130 C 

0131 URITE(40»1000) NCCNP 

0132 HRITE(40rlOOO)NAiNMfNCONS»IADDFHIRIHG»NDPSI»NPR 

0133 «RITE(40f1004) PSIFFfREVfCTfXHU.ALPHATfRBfH 

0134 URITE(40f1001)ALFRfBETRfALFZfBETZfEPS 

0135 URITE(40Fl001)((A(lFJ)Fl=lFHyi)fJ=lF2) 

0136 9RITE(40Fl002)((GAhA(lFJ)Fl=lFNy)FJ=lF2) 

0137 9RITE(40Fl00l)((XtI»J)FY(lFJ)FZ'.lFJ)Fl=lFNyi)FJ=lF2) 

0138 yRITE(40Fi001>((SX(lFJ)FSY(lFJ)»SZaFJ)Fl=lFNyi)Fj=lFNA) 

0139 STOP 

0140 C 

0141 C THIS SUBROUTINE USES THE INDUCED VELOCITIES AND 

0142 C OLD LOCATION OF EACH POINT AND CONFUTES NE« 

0143 C LOCATION OF THE POINTS ON THE BAKE USING SIHPLE 

0144 C FORWARD EULER HETHOD, 

0145 C 

0146 SO CALL NEWLOC 

0147 GO TO 20 

0148 C 

0149 C END OF THE NAIN LOOP 

0150 1000 F0RNATa2I5) 

0151 1001 F0RNAT(8F9,4) 

0152 1002 F0RHAT(7F10,6) 

0153 1003 F0RHAT<3F10.5fI10) 

0154 1004 F0RHAT(F12.2f2XfF6.3f2F12.5f3F7.2) 

0155 1005 F0RNATa8FE18.5F3E15.5) 

0156 END 


OF poOf^ 


" V 
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6“0cl-l?S3 03:00!40 WX-11 FORTRAN V3.1-23 

26-Jul-l?S3 00 : 23:31 SY3}USER:CSABERLREP3INI«hK.FOR>3 


0001 

C 

«tm SUE. inuake ttm 

0002 

C 


0003 

C 

THIS SUBROUTINE CONFUTES THE INITIAL UAKE LOCATION. 

0004 

C 


0005 


SUBROUTINE INUAKE 

0004 

C 


0007 

c 

THIS SUBROUTINE CONFUTES THE C00RDINATF.3 OF INITIAL - 

0008 

c 

WAKE AND ITS PROPERTIES . INITIAL UAKE IS A HELICAL - 

0009 

c 

— - WITH CONSTANT RADIUS, 

0010 

c 


0011 


CONNON/ ALLSUB/ ALFO f ALFHAT r C AT I C AT2 ( C A2T 1 CT f DALFO » DPSI 1 EPS » 

0012 


1 H>IADDfIAyFHIRIHG»NAfHCONS»NDPSI.NFR»NRHiNU»Nai»PI» 

0013 


2 PSI(PSIF»P3I0»RAD»RCiP.EV*SATFSAT2*SA2TfTHCAT» 

0014 


3 TFRC F THSAT F XHU F X>: F Y Y F ZZ F ALFR F BETR F ALFZ f* ETZ ; 

0015 


4 AaS0F2)FGAHA(180F2)FGANB(30)FSEGaSOf2)FUU8OF2)F 

0014 


5 VaS0F2)FU(180F2)FXaS0F2)FYaB0F2)FZ(180F2) 

0017 


C0R=.003 

0018 


IF(HIRIHG.NE.O) H=2,5 

0019 


SIE=0, 

0020 


DO 10 J=l»2 

0021 


IF(J,EQ,2) SIE=PI 

0022 


DO 10 I=1fNW1 

0023 


A(IfJ)=COR 

0024 


X(IfJ/=COS(SIE) 

0025 


Y(IfJ)=SIN(SIE) 

0024 


Z(I.J)=-(H-C0R)*(I-1)/NU1 

0027 


SIE=SI£-DPSI 

0023 


10 CONTINUE 

0029 


CALL SHOOTH 

0030 


DO 30 J=1f2 

0031 


K=NA 

0032 


IF(J.EQ.2)K=NA/2 

0033 


DO 30 1=1 fNU 

0034 


GAHA(IfJ)=GAHB(K) 

0035 


K=K-1 

0034 


IF(K.EQ.O) K=NA 

0037 


30 CONTINUE 

0033 


RETURN 

0039 


END 


ORIGINAL PAGE ® 
OF POOR QUALITY 
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«-0et-l?S3 03I01!32 WX-11 FORTHmN V3.1-23 

26-Jul-l?33 22M4I1? 3ySJUSER!CSHB£RI,REP]ADDTRH,FOR»9 


p£ 


0001 

0002 

0003 

0004 

0005 

0006 
0007 
0003 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 
0013 

0019 

0020 
0021 
0022 

0023 

0024 
0023 
0026 
0027 
0023 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 
OOw/ 
0033 


SUBROUTINE ADDP(NAHP) 


C — THIS SUBROUTINE ADDS HALF A TURN TO THE HAKE UHEN HORE- 
C ACCURACY FOR THE HAKE IS REQUIRED. THIS SUBROUTINE IS - 

C CALLED UHEN IADD=1 

C 

COHHON/ALLSUB/ALFO»ALFHAT.r/T,ChT2»CA2TfCTiDALFO»DPSI,EPS> 

1 HfIADDiIAVrHIRIHGiNAfNLC.N3fNDPSIiNPRfNRUrHU>NUhPI> 

2 P3I»P3IFfPSI0fRAD»RC>REV>SATrSAT2iSA2T»THCAT> 

3 TFRCfTHSAT»XNU»XX»YYfZZfALFF;>BETR>ALF2jBETZ> 

4 AaS0>2)iGHNA(lS0i2)>GAKB(30)»3EGa3Oi2}>U(ie0r2}f 

5 V(lSO»2)»U(lSO»2)»XaSO»2)fYa80f2)»ZaSO»2) 
COHHON/HIS/ SXUS0il3)»SY(130fl8)iSZ(lS0rl3) 
ilAH=NAHP-l 

no 10 KK=liNAH 
iiU--NWl 
HU1=HUH1 
AOWhl)=A<NU»l) 

A{NU1.2)=A(NU»2) 

GAtlA(NU»l)=BAHA(NU-l>l) 

GHNA(NUt2)=GAHA(NU-l»2) 

K«A=NH1-NAH 

HUB=NH1-HA 

X(NUM)=2,*X(NUA»2)-X(NUB»1) 

Y(H'Jhl).-2.<Y(NUA»2)-Y(NUB»l) 

Z(NUlil)=2.«(NUAf2)-Z{NUB»l; 

X(NU1»2)=^24X(NUAf1)-X(NUB»2) 

Y(Nai»2)=2.*Y(NUA»l)-Y(NHB»2) 

Z{HUl(2)=2..i:Z(HaA»l)-Z(NHB»2) 

DO 10 J=1fNA 
K=JtNAH 

IFIK.GT.NA) K=;K-NA 
SX(NUl»J)=2.»SX(NUArK)-3X(NHF.»J) 
SY(.Nai»J)-2.tSY(NHAiK)-SY(NUB>J) 
SZ(NUhJ)=2.)SZ(NUAFK)-SZ{NUBiJ) 

10 CONTINUE 




RETURN 

END 


f 


w 



6-0cl-1983 03!li:i0 VAX*U FORTRAN V3.1-23 
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0001 


SUBROUTINE VELOC 


0002 

C 

— 


0003 

C 

FOR A GIVEN WAKE LOCATION THIS SUBROUTINE CONFUTES THE 


0004 

C 

TOTAL INDUCED VELOCITY FOR ALL THE POINTS THE WAKE. 


0005 

C 



0006 


C0HN0N/ALLSUB/ALF0rALPHAT»CAT»CAT2»CA2T»CT»DALF0iDPSlFEPSf 


0007 


1 H*IADIifIAV»NIRINGfNA)NCONSrNDPSI>NPRrNP.UrNU;NWl>PI 


0008 


2 PSI»FSIF»PSI0)RAri»RCfREV»SAT»SAT2>3A2T»THCAT» 


0009 


3 . TFRCfTHSAT»XHU»XX»YYfZ2fALFfiffETfi»ALF2»BETZ» 


0010 


4 A(ie0>2)»CANA(180>2)iGANB(30)>SEG{ie0r2)fU(180»2)r 


0011 


5 V<lSO;2)»U(130,2)>X(180»2)fY(lSO»2)fZ(180»2) 


0012 


COMMON /NTURN/ NCONP 


0013 

C 



0014 

C 

TRANSITION VELOCITY 


0015 

c 



0016 


XHUX=XNUICAT/CT 


0017 


XMUZ=-XMUISAT/CT 


0018 


DO 10 I=1»NC0NS 


0019 


DO 10 J=b2 


0020 


U(I»J)=0. 


0021 


va»J)=o. 


0022 


10 wa;j)=0. 


0023 


NSTOP=NCONS4NCOMP 


0024 


IF(HST0P ,GT. NU) NST0P=NM 


0025 

c 



0026 

c 

— - START OF THE MAIN LOOP 


0027 

c 


• 

0028 


DO 28 IsNCONSfNSTOP 


0029 


13 DO 29 J=h2 


0030 

c 



0031 

c 

CHECK FOR HOVER CONDITION 


0032 

c 



0033 


IF(XNU,GT.0.00001.0R,J.NE.2) GO TO 150 


0034 


U(If2)=-U(M) 

T 

0035 


V(Ii2)=-V(M) 

n . 

0036 


U(I»2}= Udil) 


0037 


GO TO 29 


0038 


ISO CONTINUE 


0039 


XX=X(IfJ) 


0040 


YY=Y(IiJ) 


0041 


2Z=Za.J) 

M 

0042 


U(IiJ)=0. 

E 

0043 


V<IfJ)=0. 

0044 


ua>j)=o. 


0045 


IDUM=I 


0046 


JDUH=J 


0047 

c 



* 0048 

c 

INDUCED VELOCITY BY THE MIRROR IMAGE OF THE WAKE AT - 


0049 

c 

POINT ilJ) 


0050 

c 


1 

, 0051 


IF(HIRIMG.EO.O) CALL IMGUAKdDUMf JDUM) 

1 

f 

0052 

c 


0053 

c 

INDUCED VELOCITY BY THE MAIN WAKE AT POINT d»J) 


0054 

c 


r .. 

0055 


CALL WAKEdDUM»JDUM) 


0056 


IFd.GT.2*NA+4) GO TO 140 


0057 

c 
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VELOC 


6-0ct-1983 03lim0 VAX-11 FORTRAN V3.1-3 
24-JU1-19B3 03126504 SYSfUSERlCSABERI,REP]TOTVEL,FOR13 

0053 C INDUCED VELOCITY BY THE INSIDE HAKE AT POINT (I»J) 

0059 C 

0060 CALL INSHKE(IDUHrJDUN) 

0061 C 

0062 C INDUCED VELOCITY BY THE BOUND VORTICES AT POINT (I»J) 

0063 C 

0064 140 CALL BOUNDE(IDUH,JDUN) 

0065 U(I>J):U(IiJ)4XHUX 

0066 U(IiJ):H(IfJ)4XHUZ 

0067 29 CONTINUE 

0063 28 CONTINUE 

0069 C 

0070 C END OF THE HAIN LOOP 

0071 C 

0072 RETURN 

0073 END 


OF POOI^ k li 


- 74 - 


C4 


Pas 


6-0ct-l983 03'.08!02 VAHl FORTRAN V3.1-23 
26-JU1-19S3 025«;48 SVS$USER'.CSAEERI.REP3NEyL0C.F0R}6 


0001 

C 

tmt SUB. NEULQC lIMt 

0002 

c 


0003 

c 

HAVING OLD MAKE LOCATION AND VELOCITIES OF ALL THE 

0004 

c 

POINTS ON THE HAKE CGKPUTES THE NEU HAKE LOCATION. 

0005 

c 

THE METHOD OF INTEGRATION IS THE FORHARD EULER 

0006 

c 

METHOD, 

0007 

c 


0009 


SUBROUTINE NEHLOC 

0009 


COMHOrl/ALLSUB/ALFO,ALPHAT»CAT»CAT2,CA2T»CTfDALFO»DPSI>EPS» 

0010 


1 H>IADD>IAV, HIRING, NAfNCONSrNDPSIrNPRiNF:y>NU»NUliPIf 

0011 


2 PSI,PSIF,PSI0,RAD,RC,REViSAT,SAT2,SA2T,THCAT, 

0012 


3 TFRC,THSAT,XHU,XX,YY,ZZ,ALFR,BETR,ALFZ,BETZ, 

0013 


4 A(ir'),2},GAHA(180,2),GAHB<30),SEG(lo0,2),U(lS0,2>, 

0014 


5 V(130,Z},U{lSO,2),X(lSO,2},YU60,2),Z<>180i2} 

0015 


COHNON /NTURN/ NCOHP 

0016 


DATA COR/0.003/ 

0017 


HSTART=HCOHSH 

0018 


NSTOPl =NSTART + NCOHP 

0019 


IFiNSTOPl ,GT. NUl) NSTOPl =NU1 

0020 

c 


0021 

c 

COMPUTES THE NEH UAKE LOCATION 

0022 

c 


0023 


32 DO 35 J=li2 

0024 


TX1=X{NC0NS,J1 

0025 


TY1=Y(HC0NS,J) 

0026 


TZ1=Z(NC0NS,J) 

0027 


DO 33 I=NSTART, NSTOPl 

0028 


TX2=TX1 

0029 


TY2=TY1 

0030 


TZ2=TZ1 

0031 


TXl=Xa,J) 

0032 


TYl=Ya,J) 

0033 


TZ1=Z(I,J) 

0034 


X(I,J)=TX2tU(I-l,J)<DPSI»CT 

0035 


Y(I,.J)=TY2+V(M,J)IDPSI»CT 

0036 


Za,J)=TZ2+H(I-l,J)»DPSItCT 

0037 


IF(HIRIHG.NE.O) GO TO 33 

0038 

c 


0039 

c 

CHECHES IF POINT (I,J) PASSES THROUGH THE GROUND. 

0040 



0041 


320 IF(X(I,J)»SAT+Z(I,J)*CATIH.GT.A(I,J)) GO TO 33 

0042 


X(I,J)=(X(I,J)tCAT-2a,J)*SAT)ICAT-(H-A(I,J))»SAT 

0043 


Z(I,J)=-<<H-A(I,J))»CAT+(X(I,J)»CAT-Z(I,J)»SAT)»SAT) 

0044 


33 CONTINUE 

0045 


351 XJ=FL0AT(J-1)*PI 

0046 


X(1,J)=C0S(P3IFXJ) 

0047 


Y<1,J)=SIN{P3I+XJ) 

0043 


Z(1,J)=0. 

0049 


35 CONTINUE 

0050 

c 


0051 

c 

USING NUMERICAL DAMPING TO AVOID NUMERICAL INSTABILITY 

0052 

c 


0053 


CALL SMOOTH 

0054 

c 


0055 

c 

COMPUTATION OF THE NEH UAKE PROPERTIES 

0056 

c 

• 

0057 

c 
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0058 DO 38 J=li2 

0059 SEGOLD^SEGUtJ) 

0060 GAK0LD=GAHAa>J) 

0061 C0R0lD=A(lfJ) 

0062 DO 37 I=2iNy 

0063 11=1+1 

0064 Dx=xaifj)-xa»j) 

0065 DY=Ym»J)-m»J) 

0066 D2=Z(IhJ)-Z{IiJ) 

0067 SEGNEy=SQRT(DX»DX+DYtDY+DZtDZ) 

0068 CORNEU=COROLD»SQRT(SEGOLD/SEONEy) 

006 ? GANNEV=GAHOLD'ISEGOLD/SEGNEy 

0070 SEG0LD=SEG(IfJ) 

0071 C0R0LD=A(IiJ) 

0072 CAH0LD=GAHAU»J) 

0073 . SEG(I)J)=SEGNEy 

0074 An>J)=C0RNEy 

0075 GAhA(hJ)=GANNEy 

0076 37 CONTINUE 

0077 DX=X(1»J)-X(2>J) 

0078 DY=Y(1»J)-Y(2»J) 

007? D2=za»J)-2(2»J) 

0080 SEGU»J)=SQRT(DX*DX+DY»DY+DZ»DZ) 

0031 38 CONTINUE 

0032 NPS=IFIX((PSI+.05)/DPSI) 

0033 LPS=HOD(NPSiNA)+l 

0034 LFD=LP3+NA/2 

0035 IF(LPD.GT.NA) LPD=LPD-NA 

0036 GAHA(hl)=GAHB(L?S) 

0037 GAHAUi2)=GAHB(LPD) 

0088 A(M)=C0R 

0039 A(1>2)=C0R 

0090 RETURN 

0091 END 




OF POOR QUAIHY 


- 76 - 


C-0CI-19S3 03!08‘,5i Vi^X-ll FORTRAN W.1-23 
25-Jul-1993 02l525n SVSIUSER’.CSABERI.REFTSLFINO.FORjS 


0001 

0002 

0003 

0004 

0005 

0006 
0007 

C 

c 

c 

c 

c 

c 

THIS SUBROUTINE COMPUTES THE TSELF INDUCED VELOCITY 
AT THE POINTS (I;J) FOR SEGMENTS SETUEEN (I-1»J)» 
AND (HliJ). 

SUBROUTINE SLFIND(IDUH»IRM>IR»IRP»LiGAM»DU»DV»DH) 

0008 


C0MM0N/ALLSUB/ALF0(ALPHATiCAT>CAT2»CA2T»CT>DALF0»DPSI»EPS» 

000? 


1 HrIADDiIAV;MIRIM6>NArNC0NS>NDPSI>NPR>NRU>NU»NlllFpIi 

0010 


2 PSI»PSIFrFSI0»RADiRC»REV»SAT>SAT2»SA2TfTHCAT» 

0011 


3 TFRCiTHSAT»XNU>XXfYY»ZZrALFRiEETR»ALFZ>BETZf 

0012 


4 AaS0r2)>GAMA(lB0>2)iGAMB(3C)iSEG(180»2)rU(180>2)r 

0013 


5 VaS0»2);«(130»2)fX<lS0»2)fY(180»2)iZ{iS0>2) 

0014 


IFdDUM.EQ.l) GO TO 260 

0015 


Xl=X(IRMrL> 

0016 


Yl=Y(lRMrL) 

0017 


Zl=ZaRMFL) 

0013 


X2=X(IR>L) 

001? 


Y2=V(IRrL) 

0020 


Z2=Z(IR»U 

0021 


X3=X(IRPfL) 

002H 


Y3=Y(IRPfL) 

0023 


Z3=Z(IRP»L) 

0024 


DX1=X1-X2 

0025 


DX2=X2-X3 

0026 


DX3=X3-X1 

0027 


DY1=Y1-Y2 

0029 


DY2=Y2-Y3 

002? 


DY3=Y3-Y1 

0030 


DZ1=ZK2 

0031 


D22=Z2-Z3 

0032 


DZ3=Z3,-21 

0033 


AL2=DXnDXl+DYl»DYltDZl»DZl 

0034 


BL2=DX2»DX2+DY2»DY2+DZ2IDZ2 

0035 


AL~S5RT(AL2) 

0036 


BL=S0RT(P.L2) 

0037 


CL2=DX3*DX3fDY3»DY3+DZ3ID23 

0038 


CL=SQRT(CL2) 

0039 


DENR=<ALtPL-CL)»(AL+BL4CL)I(EL+CL-AL)»(AL+CL-BL) 

0040 


IF(BEHR.OT.O.OOI) go to 10 

0041 


DU=0, 

0042 


DV=0, 

0043 


DU=0, 

0044 


RETURN 

004S 


10 R2=AL2iBL2»a2/DENR 

00*1^ 


R=SQRT{R2) 

0047 


Sl5*S4:AL/R 

0048 


§2=,5IBL/R 

004? 


Cl=saRT{l,-Sl»31T 

0050 



0951 


Tiwii/a.fCi) 

0052 


204=82/(1.402) 

Ou53 


U ’-GAHArRHiL)»(AL0G<8.*R|Ti04/A(IRM?U 5f ,25) 

0054 


Ok eAHAiI/f?L))|t{AL0G(8,«*T204/AaR;L))+,25) 

0055 


Gg&*(GGlF8B2)/(4.tR) 

0056 


20 XiX=DYl«iZ2'-»Y2ID21 

0057 


XHY-DZl*DX2-rf22«Xl 


oRiGm/iL Pi':"'; a 
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0058 XHZ*I)XUDY2-DX2«Yl 

0059 XHL»S0RT(XHX«XHX+XHYIXHY+XHZIXH2) 

0060 BU»GGG»XHX/XML 

0061 DV«GGGIXHY/X«L 

0062 Dy>GGGtXHZ/XHL 

0063 CmO URITE(4SiCllll>IrJ>DUrDViDy 

0064 Cim F0RHATI28X» 'SELF INDUCED VEL0C'/2I10»3F12, 6) ,,, 

0065 RETURN ORISi«?i. S 

0066 260 DU*0* OF POOR » V 

0067 00=0. 

0068 DU=-GAHA(1>U»TFRC 

0069 C1112 URITE(4S»C1111)I>J>DU>D0>DU 

0070 RETURN 

0071 END 


6-0ct-19S3 03!0?:S7 UAX-11 FORTRAN V3.1-23 Pas 

26-Jul-1983 03106153 SYSIUSER1CSABERI.REP3SHALH.F0R15 

0001 SUBROUTINE UEUGAHiASQfDUrDU.DU) 

0002 C 

0003 C THIS SUBROUTINE STORES THE TINE HISTORY OF THE UAKE 

0004 C FOR HALF OF A REVELOTION IN SX SY S3. IT ALSO ADDS 

0005 C THE NUHERICAL DAMPING TO AVOID NUMERICAL AND UAKE 

0006 C INSTABILITY. 

0007 C 

0008 COMMQN/VELO/DXA f DYA »DZA»DX »DY f DZ »DXB.-DYB »DZB 

0009 RlS=DXAtDXA4DYAIDYA+DZAIDZA 

0010 R1=SQRT<R13) 

0011 R2S=DX»DXfDYIDY+D2tDZ 

0012 R2=SQRT(R2S) 

0013 SQ=DXBtDXBIDyBI'DYB4DZ6tDZB 

0014 3QDUM=(R1+R2)»(R1+R2)-SQ 

0015 HS0=.25»SQDUNI(SQ-{R1-?.2)I(R1-R2) )/SQ 

0016 GG=GAMI(R1tR2)/(RUR2»SQ£UM) 

0017 GGG=GG«(HSQ/(ASQ4HS0)) 

0018 161 CU=GSGI(DYAtDZB-DZAtDYB) 

0019 DV=GSG<(DZA>OXB-DXA»PZB) 

0020 DU=GGGI(DXA»DYB-DYAIDXB) 

0021 RETURN 

0022 END 


Pa 


i-0ct‘l?S3 03103130 %'Ml FCRTRCiN W4-23 
26-JUH9S3 03J33I4? SYSIUSERiCSABERI.REPJMAINUK'.FORH 


0001 

0002 

0003 

0004 

0005 
0004 

0007 

0008 
000 ? 
0010 
0011 
0012 

0013 

0014 

0015 
0014 

0017 

0018 
COl? 
0020 
0021 
0022 

0023 

0024 

0025 
0024 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 
0044 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 
0054 
0057 


C mu SUB. HAKEF 

C 

C VELOCITY INDUCED BY WAKE ITSELF. 

C 

SUBROUTINE UAKE(IiJ) 

C 

C — THIS SUBROUTINE CONFUTES THE INDUCED VELOCITY OF THE 

C — MAIN WAKE AT THE POINT (IfJ) ON THE WAKE ITSELF, 

C — (POINT KX,YY»2Z OR Xdf J)»Y(I,J)i2(I»J)) 

C 

C0HHCN/ALLSUB/ALF0»ALPHATiCAT,CAT2.CA2T»CTfDALF0»DPSI>EPSr 

1 HdAOD,IAV,hIRIHG.NA,NCONSiNDPSI»HPRfNR«»NW»NWhPIr 

2 PSI»FSIFFpSI0jF;ADrRC,REV>SAT»SAT2»SA2TdHCAT» 

3 TFfiC>THSAT,X«U,XXfYYrZZiALFRfBETRrALFZfBEfZ» 

4 A(lE0,2}fGAHA(ie0>2),CAHB(30)fSEG(180i2)iU(lS0>2)t 

5 V(lB0.2)FW(lS0.2),X{l60f2)iY(180f2)fZ(ie0F2) 
C0HN0N/KALFS/RRFXlUFYlUFZlUFX2UFY2UFZ2Uf 

1 X1fX2fX3fY1»Y2jV3fZ1fZ2fZ3 

COM/VELO/DXhfDYAfDZAfDXfDYfDZfDXBfDYBfDZB 
NH=NA/2 
TH’2,IH 
C 

C LOOP FOR THE WAKES PRODUCED BY TWO BLADES 

C 

DO 25 L=If 2 
LLsL 

DXA=XX-X(1fL) 

DYA^YY-YUfL) 

b:a=zz-z(1fL) 

RlS=riXA»DXAtIiYA»DYA4DZA*DZA 

R1=SQRT(R13) 

C 

C — LOOP FOR ALL THE PCOINTS ON THE WAKE — 

C 

DO 80 IRsIfHW 

IRP=IR+1 

hx=xx-x(irpfL) 

DY=YY-Y{IRPfL) 

DZ=ZZ-Z(IRPfL) 

R23=DX«X4DY»DY4DZIDZ 

R2=3QRT(R2S) 

C 

C — CHECK IF THE POINT (IfJ) IS END OF A SEGMENT BETWEEN 

C POINTS (IRfL) and (IR+2FL) 

IF(L.NE.J) GO TO 50 
IF(IR.LT.M.OR.IR.GT.I) GO TO 50 
IFd.EO.l.OR.IR.EQ.I-DGO TO 20 
IFdR.EQ.I) GO TO 70 
20 IRK=I-1 

IFd.EQ.l) IRH=1 
IRldRHd 

irfi=if:i+i 

IHUi1=I 

GArtHGAHAdSHFU+GAHAdRl fL) )/ 2 , 

LDUH=L 
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24-Jul-l9B3 03i33H9 SVSiUSERItSABERI.REPlHAINUK.FORH 

58 C COHPUES SELF INDUCED VELOCITY FOR THE TWO SEGHENTS 

59 C AROUND POINT (I»J). 

50 

51 CALL SLFIND(IDUH»IRH»IRl»IRPl»LDUNjGAH»DUniV>Da) 

52 10 UafJ)=U(I»J)+DU 

53 va»j)=v(bj)+Dv 

54 «{i>j)=yaij)!Du OF poo?< yy/;y«Y 

55 Cllll MRITE<45(C2221)L»IR.DU»DV»Dy 

56 C2221 F0RHAT(4X>'UAKE 'f2I8>3F12.6) 

57 GO TO 70 
53 C 

s9 C CHECK IF THE POINT (IiJ) IS VERY CLOSE TO OR VERY 

70 C FAR FROH SEGKENT BETWEEN POINTS (IRfL) AND 

n C {IR+bD 

72 C 

73 50 IF(R1.GT,TH.AND.R2.GT.TH) GO TO 70 

n C 

75 C 

76 C 

77 IF(R1.LT.0.2.0R.R2.LT.0.2)G0 TO 67 

78 68 SQ=SEGiIR.LXSEG(IR»L) 

■9 SQDUH=(RUR2>t(Rl+R2)-SQ 

:0 IF(SaDUM»LT.0,001) SQDUN*,001 

-:i HSQ=,25ISQDUNt(SQ"(Rl“R2)l(Rl“R2) )/SQ 

32 ASQ=A(IR»L)«A(IR»L) 

33 GG=GAHAUR!L)»(R1+R2)/(R1IR2PSQDUH) 

34 GGG=6GI(NSQ/(ASQPHSQ)) 

15 DXB=X(IR»L)-X(IRP.L) 

36 DYB=Y(IR.L)-Y(IRF»L) 

37 DZB=Z(IR)L)“Z(IRPfL) 

33 141 XNU1=DYAIDZB-DZA*DYB 

39 XNl!2=DZAJDXB-DXAtDZB 

^0 XNU3=DXA»DYB-DYA»DXB 

-'1 DU=XNU1»GGG 

•2 DV=XNU2«GGG 

=3 Dy=XNU3*G3G 

74 UaiJ)=U(IfJ)+DU 

•’5 Va»J)’V(IrJ)+DV 

’6 H(I»J)mU(I»J)+D« 

?7 C1112 URITEn5>C2221)L»IRiDU.DV»Da 

•S GO TC 70 

C 

.0 C POINT a»J) IS CLOSE TO THE SEGMENT, 

•'1 C 

2 67 IRH3=IR 

3 IR3=IR+1 

4 IRP3=IR3+1 

5 IF(IRP3.GT.HU1) GO TO 68 

6 LL=L 

7 C 

3 C CALL HALFST SUBROUTINE TO REDUCE THE STEP SIZE 

9 C 

0 CALL 'riALFSTtIRH3»IR3»IRP3>LL) 

1 DXA=XX-X1 
DYA=YY"Y1 
DZA=ZZ-Z1 
DX=XX-X1U 


4 
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WAKE 



0115 


DV=YY-Y1U 

0116 


DZ=ZZ"Z1U 

0117 


DXBsDX'DXA 

0118 


DYE=DY-DYA 

0119 


DZB=DZ-D7A 

0120 


GAH=GAHAaRH3»LL) 

0121 


ASQ--A{!RH3>LL>IAUR«3fLL) 

0122 


CALL UEL<GAH»ASQ>DUl»IiyijIiWl) 

0123 


DXA=DX 

0124 


DYA=DY 

0125 


D2A=DZ 

0126 


DX=XX-X2 

0127 


ItY=YY-Y2 

0123 


DZ=ZZ-Z2 

0129 


DXB=0X-DXA 

0130 


DYB=IiY-HYA 

0131 


D2B=DZ-DZA 

0132 


CALL VEL(GAH»ASQiDU2»DV2.DW2) 

0133 


U{IjJ)=Ua»J)+DUUBU2 

0134 


V(IfJ)=V(I»J)FDVUIiV2 

0135 


Ha>J)=U(I»J)+DWl+DW2 

0136 

C1113 WRITE(45rC2221) LflRrDUlrliVlr 

0137 


70 R1S=R2S 

0138 


R1=R2 

0139 


DXA=DX 

0140 


liYA^DY 

0141 


DZA=DZ 

0142 


80 CONTINUE 

0143 

C 


0144 

C 

END OF THE MAIN LOOP 

0145 

c 


0146 


25 CONTINUE 

0147 

c 


0148 

c 

END OF THE BLADE LOOP 

0149 

c 


0150 


RETURN 

0151 


END 


i-oci-iTSS 03:::::. 
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f ' "g’ 

OP POOR QUAU'sy 




w 


IliGUAK 


0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 
0073 
0079 

0030 

0031 

0032 
0083 

0034 

0035 

0036 
0087 
0038 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 
0107 
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U(I»J)=U(I»J)«U 

V{IiJ>=V(IfJ)4DV 

«(I»J)=W(I,J)+DU 

cun WRITE{45(C2221)LL»IRR»DU»DV>Dg 
C2221 F0RHAT(2X»2I7»3F12.6»' HIRING') 

GO TO 1335 

1332 IRH=IRR , . , 

IR=IRR+1 

IRP=IRH OF POOR 

IFURP.GT.Nyi) GO TO 1310 

L=LL 

C REDUCTION OF STEP SIZE IN SEGHENT LENGTH 

C 

CALL HALFST(IRK»IR»IRP»L) 

DXA=XX-(X1ICA2T-Z1»SA2T-THSAT) 

DVA=Yy-Yl 

DZA=ZZt (X 1 <SA2T+Z1 ICA2T+THCAT ) 

DX=XX-(X1U:»CA2T-Z1UISA2T-THSAT) 

DY=YY-Y1U 

DZ=ZZ+ ( X1UISA2T+Z1UKA2T+THCAT) 

DXB=DX-DXA 

DYE=DY-[iYA 

Ii2B=DZ-DZA 

GAH=-GAHA(IRH»LL) 

ASQ=A(IRHiLL)»A(IRH»LL) 

ASQ=,0001 

CALL VELIGANfASOrDUhliVlrDUl} 

DXA=DX 

DYA=DY 

DZA=DZ 

DX=XX-( X2ICA2T-Z2»3A2T-TKSAT ) 

DY=YY-Y2 

DZ=ZZH X2*SA2T+Z2»CA2Tf THCAT ) 

DXB=DX-DXA 

DYB=DY-DYA 

DZB=DZ-DZA 

CALL yEL{GAHrASQ»HU2»DV2fIig2) 

U(IfJ)=U(I»J)FDUl+DU2 

va>j)=v;iiJ)«vi+DV2 

y(iij)=u(i»j}FDyifDg2 

C1112 yRITE(45iC2221) LLfIRR»DUl»DVl»DyiiLL»IRR»DU2fDV2»Dy2 
1335 CONTINUE 
R1=R2 
XI=XJ 
YI=YJ 
ZI=ZJ 

134 CONTINUE 

135 CONTINUE 
RETWN 
END 
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0001 

0002 

0003 

0004 

0005 
0004 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 
0014 

0017 

0018 
001 ? 
0020 
0021 


C mtt SUB HIRINGF ««» 

C 

C VELOCITY INCUCEB BY HIRING (GROUND EFFECT). 

C 

SUBROUTINE IHGUAKdrJ) 

C 

C THIS SUBROUTINE CONFUTES THE INDUCED VELOCITY OF 

C THE NIRROR IHAGE OF THE THE WKE AT POINTS (I»J) 

C 

CONHON/ALLSUB/ALFO(ALPHAT>CATtCAT2»CA2TFCT»DALFO»DPSI»EPSr 

1 KrIADDdAVrNIRINGrNA>HCONS>NOPSI»NPRfNRUrNU»NUl>PIr 

2 PSI»PSIF»PSI0)RAD»RC»REV»SATiSAT2.SA2T»THCATi 

3 TFRC » THSAT f XNU »XX r YY » ZZ ; ALFR » BETR » ALFZrBETZ , 

4 A(lSOf2);GANA(lSO*2)»GANB(30)»SEG(180j2htJ(180i2)» 

5 V(lS0r2)»U{180,2)rX(180f2)fY(180f2}iZa80»2) 
CCNN0N/HALFS/RR;XlU»YlU»ZlUiX2U»Y2U»Z2U» 

1 Xl;X2»X3fYhY2fY3rZl»Z2»Z3 

COHHON/VELO/DXA » DYA > DZA »DX > DY » DZ »DXB »DYB rDZB 
C 

C LOO? FOR HAKE OF BLADE NO.l AND NO. 2 

C 


0022 130 DO 135 LL=b2 

0023 XI=X { 1 ? LL ) »CA2T-Z ( 1 » LL ) 4 SA2T-THSAT 

0024 YI=Y(1>LL) 

0025 ZI =-( Z ( 1 f LL ) 4CA2T+X ( 1 >LL ) 4SA2T+THC AT ) 

0024 RlS=(XX-XI)*(XX-XIH(YY-Ym(yY-YI)f(ZZ-ZI)KZ2-ZI) 

0027 R1=SQRT(R1S) 

0028 C 

0029 C LOOP FOR ALL THE POINTS ON THE WAKE 

0030 C 

0031 131 DO 134 IRR=bNH 

0032 IP=IF:R+1 

0033 XJ=X ( IP>LL) tCA2T-Z( IP f LL) »SA2T-THSAT 

0034 YJ=YIIPfLL) 

0035 Z J=- ( Z ( IP » LL ) ICA2T+X { IP r LL ) 4SA2T+THC AT ) 

0034 R2S={XX-XJ)»(XX-XJ)+(Y', •YJ)»(YY-YJ)+(ZZ-2J)*(ZZ-ZJ) 

0037 R2=SQRT(R2S) 

0033 IF(R1.GT.THCAT.AHD.R2.GT.THCAT) GO TO 1335 

003? IF(R1.LT.0.2,0R.R2.LT.0.2) GO TO 1332 

0040 C 

0041 C CHECK IF POINT (IiJ) IS TOO CLOSE TO OR TOO FAR FRON 

0042 C NIRROR IHAGE OF SEGNENT BETWEEN POINTS (IRRR-1,L) 

0043 C (IRRjL) 

0044 C — . • . 

0045 C 

0044 1310 SQ=SEG(IRR»LL)*£EG(IRRrLL) 

0047 SQDUH=(R1+R2)»(R1tR2)-3Q 

0043 HSQ=.25»SaDUH»(SCl-{Rl-R2)l(Rl-F;2))/SQ 

0049 ASQaAURRjLDWdRRfLL) 

0050 132 GG=-GAHAdRR>LL)»(Rl+R2)/(Rl»R2*SQDUH) 

0051 GGG=GG»:(H3Q/(ASQ+HSQ)) 

0052 133 XNUl=(YY-Ym(ZI-ZJ)-(22-ZI)l(YI-YJ) 

0053 XNU2=(ZZ-ZI)»(XI-XJ)-(XX-XI)»(ZI-ZJ) 

0054 XNU3=(XX-XI)»(YI-YJ)-{YY-YI)*(XI-XJ) 

0055 DU=XHU1*GGG 

0054 DV=XHU2IGGG 

0057 DU=XNU3*GGG 


ORIGt«i. Ff -S-: E 
OF POOR QUAUsV 


- 8 :^ 


6-0et'19S3 03!02I42 WX-11 FCRTKiiN V3.i-23 

26-Jul-19S3 C0:0S;00 3YSfUSER:CSHEERI.REP]FNDVTX.F0H-}6 


E nm SUB. BOUHCE imt 

c 

C VELOCITY INDUCED BY BOUND VORTICES. 

c 

SUBROUTINE BOUNDEdfJ} 

C 

C THIS SUBROUTINE CONFUTES THE INDUCED VELOCITY OF THE 

C BOUND VORTICES. 

C 

C0NN0H/ALLSUB/ALF0rALPHATfCATrCAT2»CA2T»CT»DALF0fDFSI>EPS» 

1 HdAIiDfIAV>HIRINGrNA>NCONSiNDPSriNPR'iNRHfNU>NUlfPI> 

2 PSI»FSIF»PSI0»RAD>RC;REVf3AT»SAT2.3A2T»THCAT» 

3 TFRC t THSAT f XNU f XX r YY »:Z » ALFR >BETR f AirZiBETZi 

4 A(180f2)i6ANA<lSCr2)>GANB(30}>SEGUBOr2>r4J(180i2)r 

5 V{130f2)»U{180.2)»Xa80i2)FY(180»2)iZa80»2) 
COHKON/VELO/DXAi DYA fDZA fDX fDY fDZfDXB iDYBfDZB 
ASQ=1./(RC*RC) 

NPS--IFIX((P3I+.05)/DPSI) 

LPS=N0IKNPS:NA)+1 

LPD--LPSINA/2 

Ir(LPD.GT.NA) LFD=LPD-NA 
IF(NIRIHG.NE.O) GO TO 10 
C 

C INDUCED VELOCITY OF THE HIRROR IMAGE OF THE BOUND VORTICES 

C 

DXA=XX+THShT 

DYA=YY 

DZA^ZZFTHCAT 

DX=XX-X(1f1)*CA2T+THSAT 

DY=YY-Y(1f1) 

DZ=ZZ+XilFl)»SA2TfTKCAT 

DXB=DX-DXA 

DYB=DY-DYA 

DZB=DZ-DZA 

GAN=-GAMB(LPS) 

CALL VEL(GAHFASQFDUlFDVlFDyi) 

DX---XX-X(1f2)»CA2TtTHSAT 

DY=YY-Y(1f2) 

DZ=Z2+X(1f2)»SA2T+THCAT 

DXB=DX-DXA 

DYB=DY-DYA 

DZB=DZ-DZA 

GAH=-GA«B(LPD) 

CALL VEL{GAHfASQiDU2fDV2fDU2) 

UaFJ)=UaiJ)+DI(l+DU2 

V(IfJ)=V(IfJ)+DVUDV2 

U(IfJ)=U(IfJ)FDH1+DW2 

cmi URITE('15FC2221)lFjFDUlFDVlFDmFlFJFDU2FDV2FDa2 
C2221 F0RHAT(5Xf 'BOUND ’f2I7f3F12.A) 

10 IF(I.NE.l) GO TO 20 
XX=.5Y.(XX+X(2 fJ)) 

YY=.5I(YYtY(2fJ)) 

Z2=.5»<22FZ(2fJ)) 


OF POOR 


INDUCED VELOCITY OF THE BOUND VORTICES 


20 DXA=XX 
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0070 
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0075 
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0060 
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DYA=YY 

DZA=Z2 

DX=XX-X(1»1) 

DY=YY-Y{lrl) 

DZ=ZZ-2(lil) 

DXB=DX-DXA 

DYB=PY-DYA 

DZB-DZ-DZA 

CAM=GAHEaPS) 

CALL VEL{GAN>ASQ;PU1;BV1>IIU1) 

DX=XX-X(lr2) 

DY=YY-Ya»2) 

DZ=ZZ“Z{1»2) 

DXB=DX-DXA 

DYB=DY-DYA 

DZB-DZ-DZA 

GAH=GAHB{LPD> 

CALL VEL(GAM»ASQrDU2fDV2»DW2) 
U(I»J)=U(IfJ)«UlFDU2 
y(I»J)=V(IfJ)+DVHDV2 
y(I»J)=«(I»J)+DWHD«2 

C1112 WRITE(45»C2221)IrJ»DUl»DVhDHljI»JfDU2»DV2»D«2 
RETURN 
END 


ORlGirJAL 

OF POOR QUALITY 
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0001 

C 


0002 

C 

VELOCITY INDUCED BY INSIDE UAKE . 

0003 

C 


0004 


SUBROUTINE IHSUKE(IfJ) 

0005 

C 


0006 

C 

THIS SUBROUTINE COMPUTES THE INDUCED VELOCITY --- 

0007 

c 

THE INSIDE WAKE AT POINTS IrJ.(Z70 OF THE BLADE 

0003 

c 

RADIUS FROM THE ROOT) 

000? 

c 


0010 


CCHM0N/ALLSUB/ALF0»ALPHAT»CATrCAT2»CA2TjCT»DALF0»DPSI»EFS» 

0011 


1 HfIADD>IAV»MIRIM6fNArNC0NS>NDPSIfNPRiNRUiNUrNWlrPI> 

0012 


2 PSIiFSIF»PSI0»RADfRC»REVfSAT»SAT2fSA2T»TKCAT> 

0013 


3 TFRC » THS AT f XMU • XX » Y Y » 22 T ALFR » BETR r ALFZ » BETZ » 

0014 


4 A(ie0r2}>GAHA(180f2)rGAMB<30)iSEGaS0r2>rtH180r2}r 

0015 


5 V(lE0»2).U(lSOi2)fX{lSO»2>»Y<13O»2).ZU8O»2) 

0016 


COMMON/ VELO/DX A » DY A f D2 A > DX » D Y , DZ ^ DXB > D YB » DZB 

0017 


NEND=3*NA/2+4 

0018 


DO 10 L=l»2 

001? 


DXA=XX-.7IX(1»L) 

0020 


DYA=YY-.7tY(lFL) 

0021 


D2A=:Z-Z(1fL) 

0022 


DO 10 IR*2»NEND 

0023 


IRM=IR-1 

0024 


DX=XX-.7»XiIR»L) 

0025 


DY=YY-,7»Y(IRiL) 

0026 


DZ=ZZ-Z(IR.'L) 

0027 


DXB=DX-DXA 

0023 


DYB=DY-DYA 

002? 


d:=dz-dza 

0030 


GHM=-»5*GAMAaRM»L) 

0031 


ASQ=A(IRHfL)*A(IRMfL)/.01 

0032 


CALL VEL(GAHFASQFDUFDyFDM) 

0033 


U(I»J)=U(I»J) 

0034 


V(I»J)=V(IfJ) 

0035 


W(IiJ)=U(I»J)fDa 

0036 

Cmi yRITE(45.C,2221)L>IRfIiUjDV»DU 

0037 

r:21 F0RHAT(2Xf' INSURE 'r2I3*3F12.6) 

0038 


DXA=DX 

003? 


DYA=UY 

0040 


DZA=DZ 

0041 


10 CONTINUE 

0042 


RETURN 

0043 


END 


mm t!1’ 

OF POOR QUALITY 


m 


w 
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0001 

0002 

0003 

0004 

0005 
OOOi 

0007 

0008 
000 ? 
0010 
0011 
0012 

0013 

0014 

0015 
OOlo 
0017 
0013 
001 ? 
0020 
0021 
0022 

0023 

0024 

0025 

002o 

0027 

0028 

0029 

0030 

0031 
0022 

0033 

0034 

0035 

0036 

0037 

0038 
003? 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 
«■ 0043 

0049 

0050 
, 0051 

0052 

0053 

0054 

0055 
0036 
0057 


C 

C 

SUBROUTINE HALFSTdRHiIRrIRPrL) 

C 

C THIS SUBROUTINE PASSES A CIRCLE THROUGH THREE POINTS 

C AND FINDS THE HALF HAY BETHEEN THREE POINTS AND 

C THE HALF 'JAY BETHEEN EACH THO POINTS. IN OTHER HORDS 

C IT REDUCES THE STEP SIZE . THIS SUBROUTINE IS CALLED 

C UHEN POINT OF INTEREST IS VERY CLOSE TO A SEGHENT OF 

C — - THE HAKE. 

C 

CUHHON/ALLSUB/ALFO » ALFHAT f CAT i CAT2 » CA2T rCT »DALFO f DPSI j EPSf 

1 HrIADD>IAVrHIRIHG>NA>NC0NS>NDP3IfNPRrNRUriiy>NyirPI> 

2 PSIfFSIF»PSI0fRAD>RC>REVfSAT»SAT2fSA2T»IHCAT» 

3 TFRCfTHSAT»XHU»XXfVY»ZZ.ALFR.BETRiALFZ»BETZ» 

4 A(180i2)»GAHA(ieOr2)»GAHB(30)»SEG(180i2)jU(lSO»2)F 

5 Va80!2)!U(180.2)»xa80»2)»Y(180»2)»Z(180»2) 
C0HH0N/HALFS/R;XlU»YlU!ZlU;X2UfY2U»Z2U. 

1 Xl»X2»X3»Yl»Y2?Y3fZl»Z2.Z3 

X1=X(IRH;L) 

Y1=Y<IRH»L) 

21=ZaRH;L) 

X2=X(IRrL) 

Y2=Y(IR»L) 

Z2=Z(IR>L) 

X3=X(IRP»L) 

Y3=Y(IRP»L) 

Z3=Z(IRPiL) 

DX1=X1-X2 
DX2=X2-X3 
DX3=X3-X1 
DY1=Y1-Y2 
DY2=Y2-Y3 
DY3=Y3-Y1 
DZ1=Z1-Z2 
. DZ2=Z2-Z3 

DZ3=Z3-Z1 

AL2=DX1IDX1+DY1JDYHDZKDZ1 

BL2=DX2IDX2+DY2IDY2fDZ2»Z2 

AL=SQRT(AL2) 

BL=SQRT(BL2) 

CL2=BX3»DX3+DY3»DY3+DZ3*DZ3 

CL-SORTICL*^) 

DENR=<AL+BL-CL)XAL+BL+CL)»(BL+CL“AL)»(AL4CL-BL) 

IF(DENR.GT.O.OOl) GO TO 10 

X1U=,5»(X1TX2) 

Y1U=.5»(Y1+Y2) 

ZlU“.5T(21+22) 

X2U=.5»(X2+X3; 

Y2U=.5t(Y2+Y3) 

Z2U=.5*(Z2FZ3> 

RETURN 

10 R2=AL2«L2*CL2/DENR 
R=SQRT(R2) 

ANUH=DX3*DX2FDY3*DY2+DZ34DZ2 
DEHl=DXl»DX2fDYl*DY2+DZl»DZ2 
BEN =DEN1»D£'N1-AL2«L2 


OF POOR QUALITY 


- 87 - 


HALFST 



U 

V 

c 


k 
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Ps 


0001 


SUBROUTINE SKOOTH 

0002 

C 


0003 

c 

THIS SUBROUTINE SHOOTHS THE LOCATION OF TIP VORTICES. 

0004 

c 


0005 


COHfjON/ALLSUB/ALFO»ALPHATiCAT»CAT2)CA2TfCT»DALFO»DPSI»EPSf 

0006 


1 H>IADIiiIAVfNIRIKQ>NAfNCONS>NDPSIrNPR)NRU>NU>NUlrPh 

0007 


2 PSI(PSIF»PSI0fRAIi>RC>REViSATfSAT2»SA2T»THCATj 

0008 


3 TFRCJHSAT»XHU»XXfYY»Z2»ALFRrBETRFALFZfBETZf 

0009 


4 A(l80f2)fGAHA(180»2)»CAHB(30)fSE8<180»2hUa80f2)F 

0010 


5 V(180»2)»W<180»2)FXa80f2)fY{180F2)»Z(18C'f2) 

0011 


COHNON/HIS/ SX{180rl8)FSY(180»18)fSZ{180»18) 

0012 


COHNON /NTURN/ NCOHP 

0013 


NSTOP = NC0N3 + NCOHP 

0014 


IF(HST0P .GT. NH) NSTOP =NW 

0015 


NST0P1= NSTOP + 1 

0016 


Nl=IFIX((PSI+.05)/DPSI) 

0017 


P8I=N1WPSI 

0018 


N2=H0n{Nl»HA) 

0019 


NAH=NA/2 

0020 


JsN2+l 

C021 


K=J+NAH 

0022 


IF(K.GT.NA) K=K-NA 

0023 


40 N3TART=NC0NS+1 

0024 


PSIE=NAH»DPSI-,1»IIPSI 

0025 


IF(PSI.LT.PSIE) GO TO 20 

0026 


IF(J,EQ.4.0R.J.EQ,10) GO TO 44 

0027 


IF{J.NE.1.AND.J.NE.7)G0 TO SO 

0028 


44 rXl=X(NSTARTFl)-SX{N3THRT»J) 

0029 


DX2=X(NSTART»2)-SX(HSTART»K) 

0030 


DY1=Y(NSTART»1)-SY(HSTART»J) 

0031 


HY2=Y(HSTARTf2)-SY(NSTARTFK) 

0032 


DZ1=Z(NSTARTf1)-EZ(NSTART»J) 

0033 


DZ2=Z(NSTART»2)-S7(NSTART»K) 

0034 


ERRl =ALFR» (DXl»DXlf DYllDYl ) + ALFZ» ( DZl »DZ1 ) 

0035 


ERR2=ALFRI(DX2»DX2FIiY2tDY2)+ALFZ»(DZ2»DZ2) 

0036 


ERR=SORT( (ERR1tERR2)/{ALFR4ALFZ) ) 

0037 


IF(ERR.GT.EPS) GO TO 50 

0038 


IF(NCONS.GE.NU) GO TO 50 

0039 


NC0N3=NC0N':'.t 

0040 


GO TO 40 

0041 


50 CONTINUE 

0042 

c 

PRINT h' 1 NSTOPl fNCOHPfNCONS 'rNSTOPlfNCOMPiNCONS 

0043 

c 

PRINT *1 ' XEfSXE 'f X(NH1f1)fSX(NH1fJ) 

0044 


DO 2 MfNCONS 

0045 


X(If1)=SX{IfJ) 

0046 


X{I,f2)=SX(IfK) 

0047 


Y(If1)=SY(IfJ) 

0048 


Y(If2)=SY(IfK) 

0049 


Z(If1)=S2(IfJ) 

0050 


Z(If2)-SZ(IfK) 

0051 


2 CONTINUE 

0052 


IF(NST0P1 ,GE, NHl) GO TO 4 

0053 


DO 3 I=NST0P1fNW1 

0034 


X(If1)=SX(IfJ) 

0055 

0056 


X(If2)=SX(IfK) 

Y(lFl)=SYaFJ) 

0057 


Y(If2)=SY(IfK) 


ORIGINAL FASC: i'y 
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00S8 


Z(Iil)=SZ{l!j) 

0059 


Z(I»2)=SZa(K) 

00^0 


3 CONTINUE 

0061 


4 CONTINUE 

0062 

C 

PRINT t»' 2 NSTOPl jNCOHPfHCONS 'fNST0PliNC0HP»NC0NS 

0063 

C 

PRINT tt ^ XErSXE X<NUlil),SX(NUhJ) 

0064 


DO 10 I=NSTART»NST0P1 

006S 


Xa»l)=ALFRIX(I»l)+B£TRISXUrJ) 

0066 


y(Iil)=ALFR4!Y{I»l)4BETR»SY{IiJ) 

0067 


X(I»2>=ALFR»X(I»2)+BErRtSXa»K) 

0068 


Y(I>2)=ALFRtY(I.2)4BETR»SY{I»K) 

0069 


Z(I»l)=ALFZ*2(Ifl)4BETZtSZ(I»J) 

0070 


Z(Ii2)=ALF2IZa»2)+BETZtSZ(IiK) 

0071 


10 CONTINUE 

0072 

C 

PRINT *>' 3 NSTOPl »NCDNP»NCONS '»NST0P1»NC0NP»NC0NS 

0073 

C 

PRINT tt ' XErSXE '» XlNMWDrSXlNWhJ) 

0074 


IFIXHU.GT.O.OOOOI.OR.HIRIHG.NE.I) GO TO 80 

0075 


Rl=l, 

0076 


DO 60 I-2>N«1 

0077 


R2=SQRT<X(Mm(Irl)+Y(I»l)IY(I»l)) 

0078 


IF(R1.GT.R2) GO TO 70 

0079 


X(M)=X(M)»R1/R2 

0080 


Y(Ifl)=Y(M)4Rl/R2 

0081 


R2=R1 

0032 


70 X(h2)=-X(I»l) 

0083 


Y!If2)=-Y(M) 

0084 


Z(I>2)= Z(I»1) 

0035 


R1=R2 

0086 


60 CONTINUE 

0087 

C 

PRINT it' 4 NSTOPl iNCOKPrNCONS SNSTOPhNCOHPrNCONS 

0088 

q 

PRINT If ' XE»SXE 'i X(N«l»l)fSX(Nai»J) 

0089 


GO TO 20 

0090 


80 CONTINUE 

0091 


CALL SH0TH2 

0092 


20 CONTINUE 

0093 


DO 30 1=1 »Ni': 

0094 


sxa»j)=i'.a:-;i. 

0095 


SX(I»K)=Xa>2) 

0096 


SY(I»J)=Y(Irl) 

0097 


SY(IrK)=Y<If2) 

0098 


SZ(IiJ)=Z(M) 

0099 


SZ(I»K)=Z(I»2) 

0100 


30 CONTINUE 

0101 

c 

PRINT tt' 5 NSTOPl iNCOHP»NCONS '>NST0P1 jNC0HPjNC0NS 

0102 

c 

PRINT tt ' XEfSXE 't X(N«1>1)»SX(N«1»J) 

0103 


RETURN 

0104 


END 


OR QU*Uuf 
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